The Hamiltonian formulation of General Relativity: myths and 

reality 

N. KiriushchevaQ and S.V. KuzmiiJ3 

Faculty of Arts and Social Science, Huron University College, 
N6G 1H3 and Department of Applied Mathematics, 
University of Western Ontario, N6A 5B7, London, Canada 



1 



Abstract 

A conventional wisdom often perpetuated in the literature states that: (i) a 3+1 decomposition 
of space-time into space and time is synonymous with the canonical treatment and this decompo- 
sition is essential for any Hamiltonian formulation of General Relativity (GR); (ii) the canonical 
treatment unavoidably breaks the symmetry between space and time in GR and the resulting al- 
gebra of constraints is not the algebra of four-dimensional diffeomorphism; (iii) according to some 
authors this algebra allows one to derive only spatial diffeomorphism or, according to others, a 
specific field-dependent and non-covariant four-dimensional diffeomorphism; (iv) the analyses of 
Dirac [Proc. Roy. Soc. A 246 (1958) 333] and of ADM [Arnowitt, Deser and Misner, in "Grav- 
itation: An Introduction to Current Research" (1962) 227] of the canonical structure of GR are 
equivalent. We provide some general reasons why these statements should be questioned. Points 
(i-iii) have been shown to be incorrect in [Kiriushcheva et al., Phys. Lett. A 372 (2008) 5101] and 
now we thoroughly re-examine all steps of the Dirac Hamiltonian formulation of GR. By direct 
calculation we show that Dirac's references to space-like surfaces are inessential and that such 
surfaces do not enter his calculations. In addition, we show that his assumption ^ofe = Oj used 
to simplify his calculation of different contributions to the secondary constraints, is unwarranted; 
yet, remarkably his total Hamiltonian is equivalent to the one computed without the assumption 
gok = 0. The secondary constraints resulting from the conservation of the primary constraints of 
Dirac are in fact different from the original constraints that Dirac called secondary (also known as 
the "Hamiltonian" and "diffeomorphism" constraints). The Dirac constraints are instead particular 
combinations of the constraints which follow from the primary constraints. Taking this difference 
into account we found, using two standard methods, the generator of the gauge transformation 
gives diffeomorphism invariance in four-dimensional space-time; and this shows that points (i-iii) 
above cannot be attributed to the Dirac Hamiltonian formulation of GR. We also demonstrate 
that ADM and Dirac formulations are related by a transformation of phase-space variables from 
the metric g^i, to lapse and shift functions and the three-metric gkm, which is not canonical. This 
proves that point (iv) is incorrect. Points (i-iii) are mere consequences of using a non-canonical 
change of variables and are not an intrinsic property of either the Hamilton-Dirac approach to 
constrained systems or Einstein's theory itself. 
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"On ne trouvera point de Figures dans cet Ouvrage. Les 
methodes que j'y expose ne demandent ni constructions, 
ni raisonnemens geometriques ou mecaniques, mais 
seulement des operations algebriques, assujeties a 
une marche reguliere et uniforme. Ceux qui aiment I 
Analyse, verront avec plaisir la Mecanique en devenir 
une nouvelle tranche, et me sauront gre d en avoir 
etendu ainsi le domaine. " 



J. L. Lagrange, "Mecanique Analytique" (1788) 

The reader will find no figures in this work. The 
methods which I set forth do not require either con- 
structions or geometrical or mechanical reasonings, 
but merely algebraic operations subjected to a regular 
and uniform rule of procedure. Those who are fond 
of Mathematical Analysis will observe with pleasure 
Mechanics becoming one of its new branches and they 
will be grateful to me for having thus extended its 
domain. 

I. INTRODUCTION 



We begin our paper with words written more than two centuries ago by Lagrange in 
the preface to the first edition of the "Mecanique Analytique" [l| because they express 
our standpoint in analyzing of the Hamiltonian formulation of General Relativity (GR). 
The results previously obtained by others are reconsidered and classified as either "myth" 
or "reality" depending on whether they were obtained by what Lagrange called a regular 
and uniform rule of procedure, or by geometrical or some other reasonings. The results 
and conclusions constructed using such reasonings must be checked by explicit calculation; 
without which they are meaningless and could be misleading, contradicting the rules of 
procedure and the essential properties of GR. 

Originating more than half a century ago, the Hamiltonian formulation of GR is not a 
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new subject. It began with advances in the Hamiltonian formulation of singular Lagrangians 
due to the pioneering work of Dirac on generalized (constrained) Hamiltonian dynamics . 

We restrict our discussion to the original Einstein metric formulation of GR. The first- 
order, metric-affine, form 31] will be just briefly touched; but the analysis presented here can 
and must be extended to a metric-affine form and to other formulations. 

In chronological order (which is also ranked inversely in popularity) the Hamiltonian 
formulation of GR was considered by Pirani, Schild, and Skinner (PSS) 0, Dirac fl, and 
Arnowitt, Deser, and Misner (ADM) Q| and references therein. The relationship among 
these formulations has not been analyzed; and some authors have adopted to using the 
name "Dirac-ADM" or refer to Dirac when actually working with the ADM Hamiltonian. 
This presumes equivalence of the Dirac and ADM formulations. These two, as we will 
demonstrate, are not equivalent. 

The Dirac conjecture j?], that knowing all the first-class constraints is sufficient to deduce 



y, 3 



and became 



the gauge transformations, was made only after the appearance of 
a well defined procedure only later The application of such a procedure to field 

theories was considered for the first time by Castellani 101 (for alternative approaches see 

~i n n 

111 . |12| . Il3l|). Deriving the gauge invariance of GR from the complete set of the first-class 
constraints should also be viewed as a crucial consistency condition that must be met by any 
Hamiltonian formulation of the theory; yet, this requirement did not attract much attention 



and it is not discussed in textbooks on GR, where a Hamiltonian formulation is presented 
15|). In books on constraint dynamics [l6, l3, Q, even if such a procedure is 



liu 11 la 
3.g. Q, 



discussed 
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it is not applied to the Hamiltonian formulation of GR. Recently this question 
was again brought to light by Mukherjee and Saha [lol] who applied the method of 12] to the 
ADM Hamiltonian with the sole emphasis on presenting the method of deriving the gauge 
invariance, not on the results themselves. In |l9| there appears a first complete derivation 
of the gauge transformations from the constraint structure of the ADM Hamiltonian. The 
expected transformation of the metric tensor is 20 ] 



where is the gauge parameter and the semicolon ";" signifies the covariant derivative. In 
the literature on the Hamiltonian formulation of GR, the word "diffeomorphism" is often 
used as equivalent to the transformation ([!]), which is similar to gauge transformations in 
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ordinary field theories. This meaning is employed in our article.^ The expected invariance 
([1]) does not follow from the constraint structure of ADM Hamiltonian and a field- dependent 
and non-covariant redefinition of gauge parameters is needed^ to present the transformations 



of 



19| in the form of ([T]), i.e. 



Ofc 



^ — \ 9 ) ^ADMi S — ^ADM + ^00 V 9 ) ^ADM- 

The field-dependent redefinition of gauge parameters ([2]) goes back to work of Bergmann 
and Komar [21] where it was presented for the first time. The same redefinition of gauge 
parameters (121), but in a less transparent form, was obtained for the ADM Hamiltonian 
by Castellani [lO^ for the transformation of the Qq^ components of the metric tensor to 
illustrate his procedure for the construction of the gauge generators. This redefinition of 
gauge parameters was also discussed from different points of view in [22|, |23|, l2J, |25|, the most 



recent derivation is in [19| . A common feature of these different approaches is that they only 



consider the ADM Hamiltonian. According to the conclusion of 2l|, the transformation 



|T]) and the one with parameters that depend on the fields ([2]) are distinct. In [23[ this 
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transformation is called the "specific metric-dependent diffeomorphism" . The authors of 
have a brief and ambiguous conclusion about ([2]): "[it will] lead to the equivalence^ between 
the diffeomorphism and gauge transformations" and, at the same time, "demonstrate the 



unity of the different symmetries involved" ; t 



Soon after appearance of 



lese are contradictory statements. 



191 ] ■ Samanta [261] posed the question "whether it is possible to 



describe the diffeomorphism symmetries without recourse to the ADM decomposition" . To 
answer this question, he derived the transformation ([T]) starting from the Einstein-Hilbert 
(EH) Lagrangian (not the ADM Lagrangian) and applying the Lagrangian method for recov- 



ering gauge symmetries based on the use of certain gauge identities that appear in 



[ill- 



It is 



important that ([T]) follows exactly from this procedure without the need of field-dependent 
and non-covariant redefinition of the gauge parameters, which would be necessary in 
where the ADM Hamiltonian is used. The question of the equivalence of ([T]) and ([2]) does 



^ In mathematical literature the term diffeomorphism refers to a mapping from one manifold to another 

which is diffcrcntiable, one-to-one, onto, with a differentiable inverse. 
^ More detail on the derivation of ^ is given in the last Section where application of Castellani's procedure 

to the ADM Hamiltonian is reexamined {e^Df^j and e^jjj^/ are gauge parameters of ADM formulation). 
^ Here and everywhere in this article the Italic in quotations is ours. 
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not even arise in the approach of 26|]. In 26| the diffeomorphism transformations were also 
derived by applying the same method to the first-order, afiine-metric, formulation [3] of GR. 



The conclusion of 



261 ] that "the ADM splitting, which is essential for discussing diffeomor- 



phism symmetries, is bypassed" contradicts the obtained result. Firstly, any feature that 
is ^^essential" cannot be ^^bypassed" . Secondly, the transformations derived from the ADM 
Hamiltonian in |19l] are not those of 26|. It is not a "bypass" because the "destination" of 
having the invariance of ([1]) is changed. 

The conclusion about the results of [3| and 26] should be that the ADM decomposition 



is inessential and incorrect because it does not lead to diffeomorphism invariance. This 
discrepancy between these two recent results vindicates Hawking's old statement 27| "the 
split into three spatial dimensions and one time dimension seems to be contrary to the whole 
spirit of relativity", the more recent statements of Pons [24^: "Being non-intrinsic, the 3+1 
decomposition is somewhat at odds with a generally covariant formalism, and difficulties 



arise for this reason" , and Rovelli [28j : "The very foundation of general covariant physics is 
the idea that the notion of a simultaneity surface all over the universe is devoid of physical 
meaning" . 

There is another statement in 26|] that can also be found in many places "it is well known 



that this decomposition plays a central role in all Hamiltonian formulations of general rela- 
tivity" . This sentence combined with Hawking's "spiritual" statement forces one to conclude 
that the Hamiltonian formulation by itself contradicts the spirit of GR. This resonates with 

n 

Pullin's conclusion |29| that "Unfortunately, the canonical treatment breaks the symmetry 
between space and time in general relativity and the resulting algebra of constraints is not 
the algebra of four diffeomorphism" . We will show in this paper that the canonical formal- 
ism is in fact consistent with the diffeomorphism ([1]) when the Dirac constraint formalism 
is applied consistently and that the discrepancies between the ADM formalism and ([1]) can 
be explained. 

The difference of the results 19|] and 261] which were obtained by different methods 
also implies the non-equivalence of the Lagrangian and Hamiltonian formulations. In all 
field theories (e.g.. Maxwell or Yang-Mills) the Hamiltonian and Lagrangian formulations 
give the same result for gauge invariance, so for GR to differ seems unnatural. Could this 
be a peculiar property of GR? Is GR a theory in which the Hamiltonian and Lagrangian 
formulations lead to different results or was a "rule of procedure" broken somewhere? 



7 



Recently, in collaboration with Racknor and Valluri [30|, we demonstrated that, by fol- 
lowm^hem^t natural first attempt of PSS fl arrd by applying the rules of procedure 
2|, |7|, llO, lla, ll8|, the Hamiltonian formulation of GR (without any modifications of the 



action or change of variables) leads to consistent results. The gauge transformation of the 
metric tensor was derived using the method of [l^ and, without any field- dependent redefi- 



nitions of gauge parameters, it gives exactly the same result as the 



261], as it should. In the Hamiltonian formulation of GR given in 



jagrangian approach of 



301 1 the algebra of con- 



straints is the algebra of "/owr diffeomorphism" , in contradiction to the general conclusion 
of 291] which was based on the particular, ADM, formulation. 

The procedure of passing to a Hamiltonian formulation in field theories based on the 
separation of the space and time components of the fields and their derivatives (defined on 
the whole space-time, not on some hypersurface) is not equivalent to separation of space-time 
into space and time. For example, by rewriting the Einstein equations in components (as was 
done before Einstein introduced his condensed notation), we do not abandon covariance even 
if it is not manifest. In addition, such explicit separation of the space and time components 
and the derivatives of the fields does not affect space-time itself and is not to be associated 
with any 3+1 decomposition, slicing, splitting, foliation, etc. of space-time. The final result 
for the gauge transformation of the fields can be presented in covariant form when using the 
Hamiltonian formulation of ordinary field theories (e.g., Yang- Mills, Maxwell), as well as in 

grH- 

In any field theory, after rewriting its Lagrangian in components, the Hamiltonian 
formulation for singular Lagrangians follows a well defined procedure. Such a procedure 
is based on consequent calculations of the Poisson brackets (PB) of constraints with the 
Hamiltonian using the fundamental PBs of independent fields. In the case of field theories 
they are 



{g(x°,x),p(y°,y)}^„^^„ = 5(x-y). (3) 

This is a local relation that does not rely on any extended objects or surfaces. Again, 
as with separation into components, this locally defined canonical PB does not affect space- 
time and is not related to space-like surface or any other hypersurface because ([3]) is zero 
for x 7^ y in a whole space-time and there is no information in (l3l) that, using mathematical 
language, can allow one to classify two separate points as points on a particular space-like 
surface or on any surface. The canonical procedure does not itself lead to the appearance of 

8 



any hyper surf aces; in [30| there are no^references to such surfaces and the result is consistent 
with the Lagrangian formulation of [2(3]. Such surfaces are either a phantom of interpretation 
or canonical procedure was abandoned by their introduction. 

The discussion of an interpretational approach is not on the main road of our analysis 
of the Hamiltonian formulations of GR. However, the routes of such an approach'^ are quite 
interesting: one starting from the basic equations of the ADM formulation, according to 



3l| , "would like to understand intuitively their geometrical and physical meaning and derive 
them from some first principles rather than by a formal rearrangement of Einstein's law" . 
By taking this approach, a formal rearrangement (which is a "rule of procedure" ) is replaced 
by some sort of intuitive understanding. As a result, a new language is created which "is 
much closer to the language of quantum dynamics than the original language of Einstein's 
law ever was" fsi)] . This language allows one "to recover the old comforts of a Hamiltonian- 
like scheme: a system of hypersurfaces stacked in a well defined way in space-time, with the 
system of dynamic variables distributed over these hypersurfaces and developing uniquely 
from one hypersurface to another" |32]. Such an interpretation, although 'reasonable' from 
the point of view of classical Laplacian determinism, is hard to justify from the standpoint 



of GR 33]. In GR, an entire spatial slice can only be seen by an observer in the infinite 
future [3J| and an observer at any point on a space-like surface does not have access to 
information about the rest of the surface (this is reflected in the local nature of ([3]) in field 
theories). It would be non-physical to build any formalism by basing it on the development 
in time of data that can be available only in the infinite future and trying to fit GR into 
a scheme of classical determinism and nonrelativistic Quantum Mechanics with its notion 
of a wave function defined on a space-like slice. The condition that a space-like surface 
remains space-like obviously imposes restrictions on possible coordinate transformations, 
thereby destroying four-dimensional symmetry, and, according to Hawking, "it restricts 
the topology of space-time to be the product of the real line with some three-dimensional 
manifold, whereas one would expect that quantum gravity would allow all possible topologies 



of space-time including those which are not product" ^]. T 



slicing of space-time, must be lifted at the quantum level 



lis restriction, imposed by the 



351 ]: but, from our point of view. 



^ We have to confess that we found hard to understand approaches which are not analytical and, to avoid 
any misinterpretations, we will merely quote their advocates. A reader interested in this approach can 
find more details in the articles we cite. 
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avoiding it at the outset seems to be the most natural cure for this problem. 

The usual interpretation of the ADM variables, constraints, and Hamiltonian obviously 
contradicts the spirit of relativity. With restrictions on coordinate transformations which 
are imposed by such an interpretation it is quite natural to expect something different from 
a diffeomorphism transformation, as was found in 10|, Il9 |. 

Any interpretation, whether or not it contradicts the spirit of GR, cannot provide a 
sufficiently strong argument to prove or disprove some particular result or theory, because 
an arbitrary interpretation cannot change or affect the result of formal rearrangements. The 
transformation different from diffeornorphism that follows from the ADM Hamiltonian is 
the result of a definite procedure 10|, ll9l] and is based on calculations performed with their 
variables and their algebra of constraints. From the beginning we will not use the language of 
3+1 dimensions, so as to avoid the necessity of getting ourselves "out of space and back into 



space-time" [36!] at the end of the calculations. In any case, it would hkely be impossible to 
do so after we have gone beyond the point of no return on such a road. We must reexamine 
the derivation of ADM Hamiltonian right from the start. 

It is difficult to compare the results of [3^ directly with those of ADM because some 
additional modifications of the original GR Lagrangian were performed by ADM and it is 
not easy to trace them according to the "rules of procedure" . We will start with the work of 
Dirac 5|, where all modifications and assumptions are explicitly stated making it possible 
for them to be checked and analyzed. In addition, Dirac's canonical variables are compo- 
nents of the metric tensor which are the same as those used in uifll where diffeomorphism 



invariance was derived directly from the Hamiltonian and constraints. Moreover, in 



QQ and fl. 



37| two 



Hamiltonian formulations, based on the linearized Lagrangians of [4, |30| and [5| , were con- 
sidered. Despite there being different expressions for the primary and secondary constraints, 
these two formulations have the same algebra of PBs among the constraints, and with the 
Hamiltonian, therefore, they have the same gauge invariance. This is exactly what one can 
expect in the case of full GR, provided one makes no deviation from canonical procedure. 
In analyzing the ADM formulation we will follow a different path. We will not start from 
the GR Lagrangian, but instead compare the final results of Dirac and ADM and try to 
determine what deviation from the canonical procedure lead to the transformations found 



m 
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19j which are distinct from those of ([T]) 



In the next Section we shall thoroughly reexamine the Dirac derivation of the GR Hamil- 
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tonian |5| with emphasis on the effect of his modifications of the action and of the other 
simphfying assumptions he makes. In particular, we will investigate whether space-like sur- 
faces actually play any role in his derivation, or if they just serve as an illustration which 
can be completely disregarded from the standpoint of the canonical procedure, as in [30]. In 
Section III, using Castellani's procedure and the results of Section II, we derive the trans- 
formations of the metric tensor. The result is the same as those found in 26|] and 30 1. 



The same result is obtained by application of the method used in |l9:] to Dirac's Hamil- 
tonian, which illustrates the equivalence of these two methods. Some peculiarities of such 
methods, that cannot be seen in ordinary field theories, are briefly discussed and related 
to the peculiarities of diffeomorphism invariance as it compares to the gauge invariance in 
ordinary theories. Finally, we consider the ADM Hamiltonian formulation of GR. In the 
last Section IV we demonstrate that the ADM formulation follows from Dirac's by a change 
of variables. The canonicity of this change of variables (the ADM lapse and shift functions) 
is analyzed. Based on this analysis, the general and more restrictive criteria for a canonical 
transformation in the case of singular gauge invariant theories are discussed. 

II. ANALYSIS OF DIRAC DERIVATION 



In [30l] the GR Hamiltonian, constraints, closure of the Dirac procedure, and the dif- 
feomorphism transformation of the metric tensor were derived without any reference to 
space-like surfaces, the use of any 3+1 decomposition of space-time, or slicing, splitting, 
foliation, etc., as well as without modifications of the Lagrangian or the introduction of any 
new variables. (The canonical variables of j3o| are components of the metric tensor.) Dirac, 



when considering the Hamiltonian formulation of GR in [5|, also used the metric tensor as a 
canonical variable; but he made frequent references to space-like surfaces. If such surfaces, 
which according to Hawking contradict the whole spirit of General Relativity, are the part of 
Dirac's calculations, then one has to expect transformations different from diffeomorphism 
and similar to the one found in 19|] from the ADM Hamiltonian. Our main interest is to 
find out, by following all the steps of Dirac's derivation of the Hamiltonian, the place where 
(if anywhere) space-like surfaces enter his derivation or where (if anywhere) his approach 
deviates from a regular and uniform rule of canonical procedure. If there is no deviation, 

n 

one should then obtain the diffeomorphism invariance ([I]), the same as found in [30|. This 
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would resemble what happens in linearized GR, as discussed in l37l . 

n 

In [5[ , Dirac started the Hamiltonian formulation from the "gamma-gamma" part of the 
Einstein-Hilbert (EH) Lagrangian (Eq. (D8))^ (e.g., see Q) 



where 



-99 



f-pp -per _ -per -pp \ 
V pa U(J ) 



-99pu,p9aP,aB^'''^P" 



^pupapa ^ (^gf^^g'^P _ g^''^g^P) gP^ + 2 (^'^^^"^ - ^''"^'^^) 9' 



The same Lagrangian was used in J] and 301]. This is a Lagrangian of a local field theory 



(4) 



(5) 



in four(or any) -dimensional space-time, and space-like surfaces or any other hypersurfaces 
are not intrinsic to such a formulation. 

The primary constraints (the 0-equations of [5]) that follow from (jl]) are 



6L 



G 



S9pO,i 



0, 



(6) 



where p'^^ are momenta conjugate to g^^. The exact form of 0^° can be found in [j, |3] 
(Greek subscripts run from to (i — 1 and Latin ones from 1 to (i — 1 where d is the dimension 
of space-time). 

In addition to eliminating the second order derivatives of the metric tensor present in the 
Ricci scalar in passing from the EH Lagrangian to its gamma-gamma part (jll) so that 



Leh = V^R = Lg + d,V^, 



(7) 



Dirac made an additional change to the Lagrangian in order to eliminate the second term 
in ([6]). The modified Lagrangian is obtained by adding two total derivatives which are 
non-covariant (Eq. (D15)) 



L* = Lg + 



(8) 



We will refer on Dirac equations quite often and use the convention, Eq. (D#), to mean equation ^ ft'om 
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This change does not affect the equations of motion, but leads to simple primary constraints 
(Eq. (D14)) 



It was shown in [37j, that the linearized version of the modified dH]) and unmodified 
Lagrangians (jl]), despite leading to different expressions for the constraints and the Hamil- 
tonian, result in the same constraint structure, the same number of first-class constraints, 
and the same gauge invariance, which is the linearized version of_diffeomorphism. This is 
what one can also expect in the case of full GR. According to [5], the simplification iQ 
"can be achieved only at the expense of abandoning four-dimensional symmetry" which is 
obviously correct for this modification of the Lagrangian ([H]); yet Dirac's further conclusion 
that "four-dimensional symmetry is not a fundamental property of the physical world" is 
too strong and has to be clarified. Of course, four- dimensional symmetry of the Lagrangian 
is destroyed by the modification ([8]); but this change does not affect the equations of motion, 
which are the same as the Einstein equations. Consequently, for the equations of motion, 
not only four-dimensional symmetry is preserved, but also general covariance.^ If four- 
dimensional symmetry is preserved in the equations of motion, which are invariant under 
general coordinate transformations, then diffeomorphism should be recovered in the course 
of the Hamiltonian analysis, as in 30(]. 

The new Lagrangian L* differs from the original one (jl]) only for terms linear in the time 
derivatives of a metric (i.e. 'velocities'), the parts responsible for the simplification of the 
primary constraints. We then have 



r = Lg(2)+L*(1) + Lg(0), (10) 

where the numbers in brackets indicate the order in velocities (for the Hamiltonian and 
constraints it will indicate the order in momenta). The exact form of L* (1) is given by Eq. 
(D18). 

This Lagrangian is used to pass to the Hamiltonian 

^ The term "four-dimensional" symmetry used by Dirac probably reflects the fact that the gamma-gamma 
part of the Lagrangian, quadratic in first order derivatives, is not generally covariant after the elimination 
of terms with second order derivatives in the full EH Lagrangian ([7]). 
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H = g^f3,oP''^-L*. (11) 

With the modification of (IH]) the part of the Lagrangian Lq (2) + L* (1), as was shown by 
Dirac, can be written as 



1 



Lg (2) + L* (1) = Lx (0) - ^ (12) 



where is the Christoffel symbol 

= ^9'"' {9au,/3 + 9l3u,a " 9a(3,u) (13) 

and 

^rsab ^ ^rs^ab _ ^ra^sb ^^^^ 

with 

= s'' - '-Jr- (15) 

Note, that in the second order formulation, F^^, E""^^" , and are just short notations and 
none of them denote a new and/or independent variable. 

Some comments about f|T2l) are in order. The careful reader will definitely wonder how 
the parts of the Lagrangian which are quadratic and linear in velocities can have contribu- 
tions without velocities, Lx (0); the direct calculation of Lq (2) + L* (1) does not have such 
contributions (see Dirac's unnumbered equation preceding (D19)) 

^V^E™^^ [9rsfl9abfl9''' + '^9rsfl9ab,.9''' ' ^9rsfl9aP,b9^''] • (16) 

Dirac completed this square, leading to the compact form of flT2|) . Working with f|T6l) instead 
of (fT2l) . will of course not change the results and actually has no calculational advantage. 
However, we keep f fT2|) so as to compare our calculations with those of Dirac. 

The Lx (0) in (IT^ (explicitly given by (12^ ) is independent of the velocities. The only 
part of ( ITTi) that has dependence on ^f^^^o is 



9rs,of' + V^9—,E^"''ns^lb- (17) 
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Performing the variation ^ , we obtain (see (D18-D21)) 



(18) 



Equation ([T8|l is easy to solve for gah,Q due to the invertabihty of E 



rsab 



rprsab j 



(19) 



where the inverse to E"^^""^ in any space-time dimension d (except c? = 2) is 



'■abmn 



d-2 



gabgmn gamgbn 



(20) 



This result gives 



(21) 



After substitution of fl?I]) into f[T7|) (note that flTS]) can be solved for F^j, thus making the 
calculations shorter) we obtain the total Hamiltonian 

Ht = goofip'' + 2gok,op"' + Hg, (22) 
where Hq (the canonical part of the Hamiltonian) is given by Eqs. (D33, D34)) as, 



Hr 



f'gr,,. - 2 {p'^'gr.) ^ - Lx (0) - Lg (0) , (23) 



with Lx (0) (Eq. (D19)) and Lg (0) (Eq. (D8)): 



(0) = -V^W^?"^^^''"^*- 



(25) 



Note that the second term of (125]) . the part linear in the momenta, arises only after some 
rearrangement. The direct substitution of grs,o into grs,oP^^ (the only part of (ITT!) that leads 
to terms linear in the momenta) gives 
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2p 9m0,n ~l~ 00^ {^9mk,n 9mn,k) 



which after integration by parts and using |oo = —9ome leads to 



(26) 



29: 



mO 



ran , mk van ■ 

",n ~r e J ymk,n 2^ymn,k 



(27) 



The first term of ( 1271) can be written in the form given by Dirac (D41) 



29 



mO 



P. 



mn _|_ mk mn 



9mk,n i29mn,k 



9m0& iLv, 



(28) 



with 



■H,=p''9rs,v-'2{p'-''9rv),s- (29) 

We note that in obtaining the expression for the Hamiltonian 0231) . all direct calculations 
with the initially modified Lagrangian (fT2|) were performed by Dirac without any reference 
to space-like surfaces or any additional restrictions or assumptions. 

The next step in the canonical procedure is to find the time development of the primary 
constraints and see if there are any secondary constraints (or x-equations in Dirac's termi- 
nology). PBs among the primary constraints are obviously zero, = 0. The PBs 
of the primary constraints ^ with the total Hamiltonian fl2^ are 
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S9i 



-Hg = X 



Oct 



(30) 



Oct 



where we keep Dirac's convention for the fundamental PB (Eq. (Dll)), 



(31) 



[p-P (x) , g,, ix') } = \ (5^5,^ + 5l5:) 6s (x - x') . 

According to Dirac, "the second term of (D33) [Hg (0) = -Lx (0) - Lg (0) in our (|23l)] 
is very complicated and a great deal of labour would be needed to calculate it directly" 
and instead of performing the variation j^Hg (0), he uses some arguments (see (D23- 
D27)) related to the displacements of surfaces of constant time, and thus he infers that the 
Hamiltonian "must be of the form" (see (D28)) 
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where Tii and Tig are independent of the (70^. Dirac's arguments are very general and 
independent of the particular form of the Lagrangian, i.e. they have no connection with 



his initial modifications of Lq leading to L* . And, even in the linearized case 37|, without 



the=e n,odfficatio,., the .eco„da.y constraints have a dependence on s„,; this dependence 
also happens in full GR [30|. In any case, the explicit form of the constraints cannot be 
found using such arguments and explicit calculations are needed; one has to use a well 
defined rule of procedure to find them, i.e. we must calculate j^Hg (0) . Dirac performed 
these calculations using an additional simplifying assumption (see below) and this result has 
to be analyzed and compared to what follows from direct calculations. 

According to Dirac, there are no contributions from Hq (0) to a vector constraint 
(7i^ = e^'^Tis in Dirac's notation) which presumably comes from the time development of 
the corresponding primary constraint 0''° = p^^ (1301) . Furthermore, Ti^, which comes from 
the time development of the primary constraint, 0°°, can be calculated with the additional 
simplifying assumption g^o = 0, which gives (Eq. (D36)): 

= 0, = g'' = — . (32) 

fi'oo 

As a result, all of Lx (0), along with the biggest part of Lq (0), is dropped from his cal- 
culations. According to Dirac [sj], the equation for T-Cl "must hold also when gro does not 
vanish" . It is important to check this assumption by direct calculation because if the result 
of j^Hg (0) is the same as that of Dirac's, then the simplifying assumption of fl32|) . along 
with any references to surfaces of constant time, has nothing to do with his final result. 
In such a case, Hawking's criticism of formulations based on the introduction of space-like 
surfaces, which is in contradiction with the whole spirit of General Relativity and restricts 
topology of space-time [27|, cannot be applied to the Dirac analysis of GR. This also means 
that the transformations ([1]) should be derivable in the Dirac Hamiltonian formulation, as 



26| or for the Hamiltonian formulation obtained in 



was done in the Lagrangian formulation 

If the results following from the assumption of (l32l) are different from those where the 
assumption is not made, then we cannot use ( l32i) as an extra condition in the midst of the 
calculations and we have to go back to the original Lagrangian to introduce this condition 
from the outset. This is the rule followed in ordinary constraint dynamics; all imposed con- 
straints must be solved at the Lagrangian level, or added to the Lagrangian using Lagrange 
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multipliers, before performing a variation and/or considering the Hamiltonian formulation. 
For example, when Chandrasekhar considers the Hamiltonian for Schwarzschild space- 



time he, first of all, writes t 
Hamiltonian formulation 



le Lagrangian using this metric and only then passes to the 
391 ] ■ Similarly, the condition fl5^ corresponds to a particular 
coordinate system, one which is static 20, |38||; and, of course, the momenta p'^'^, which are 



conjugate to the eliminated variables gok cannot appear in such a formulation. Note that 
the initial modification of the Lagrangian ([8]) is irrelevant in a static coordinate system as 
the last two terms in ([8]) are zero when g'^^ is zero. For field theories, especially generally 
covariant ones, there is an additional restriction: the unambiguous canonical formulation 
must be performed without explicit reference to ambient space-time by making an a priori 



40|, i.e. without 



choice of a particular coordinate system or subclass of coordinate systems 
destroying the main feature of a theory from the beginning. 

To find out whether or not Dirac's formulation is correct or any reference to surfaces 
of constant time and the simplifications of [or (D36) of [5]] are relevant to his actual 
results, we perform a "great deal of labour" to find the functional derivatives separately 
for each contribution of Hq (0) = —Lx (0) — Lq (0) and to compare the results with those 
obtained by Dirac. 

For Lg (0) in ([25]), we find that 

Xg (0) = {/^ -Lg (0)} = -^V^gapM {g^'^E^^^' - g^'E^^^'^ - g^PE'^^^) 

+ \^gg,u,kgap,t [C"^'^'"^* (eee) + C^^'^^' (ee)] , (33) 

where the C's are combinations of terms of different order in e"'^ (note that the terms of 
first and zero orders in e"^ cancel) 



and 



00/30 
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When (7 = 0, the result fl33l) . is considerably simplified (this is because or E^"'^^ equal 
zero if at least one index is zero): 



Xg (0) = -lV^gaPM9''E'^('''+^^g,,,kgap,t 



V 2 



(36) 



According to Dirac, this (0) is the only source of contributions to the scalar constraint 
and he constructed it using the simplifying assumption of ( l32l) and later concluded that 
it "must hold also when qq^. does not vanish". Let us check this assertion by explicitly 
separating all space and time indices in (l36l) 



+ lV^gmn,kgp,,t (^-^E'^^P'^g'^g'^ + E'^^-^^fO^^-o + i^;-*'=P^90^no ^ 2Ef^'=-^«o^*o^ . (37) 

Some terms in fl37|) have explicit dependence on the space-time components of the metric 
tensor and these components will disappear only if condition (1321) is imposed. For (0) 
there are even more such components. Even with condition ( l32l) . the result is not zero and 
this part of the Hamiltonian, Lq (0), contributes to the vector constraint. 

Now let us find contributions coming from the second part, Lx (0). After a rearrangement 
of the terms given in into a form which is more suitable for calculation, we obtain 



Lx (0) 



2 g' 



00 



rprsab^ ^ ^fO^tiO „ „ /' rprsab , rprsba\ ^uO ^00 

-E grs,ugab,vg g - gai3,bgrs,u [J^ + ^ ) g g 
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For the 



'^'^^^"^ part we calculate 



(38) 



X'x{^) = {p"',-Lx (0)} 



2 v^*^: {-grs^ubE^^-'g'"'' + gr.,si>E'"''g'''') + (eee) + C7 (ee) + C/, (ee) . (39) 

The variation obviously produces contributions which are only third and second 

order in e"^ as in (1331) . For terms of third order we find 



(eee) = — ^ 



-ggf^u,kgaf3,t 



X 



and in second order we have two contributions: the first proportional to g^'^ 



1 ^ 



^grs,ugat,vE^""'g'''g'" - gap,,g^' {E^'^' + E^''"") (^?.s.«^?"° - grc^sg"-') 

(4 



and the second with an index a on E"^^"^ 



Ch (ee) 



-ggtMu,k- 



lOO 



' « TPrsak I rprska\ AO /30 / rpcrkrt , rpcTktr\ 

:grs,t [E + E ) g - grp,tg [e + e ) 



(42) 



Note, that we cannot present the part quadratic in e"^ f l4Tl l42l) in a compact form, where 
terms with derivatives are a common factor, because of the mixture of four and three indices, 
which is the result of the original noncovariant modification (IHl) of the Lagrangian. When 
performing these calculations we have to consider all possible combinations separately. 

It is not difficult to confirm that Xx (0) zero, even with assumption (l32l) . there are 

contributions to both the constraints Xx (0) ^^'^ (0)- Consequently, Dirac's conjecture, 
if made separately for Lx{0) and Lg{0), is not correct; but, when both parts are combined, 
the contribution of zeroth order to the secondary constraint is greatly simplified 
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x'" (0) = x'S (0) + X'x (0) 



(43) 



The (O)-part is the same with or without condition fl32l) and ~)^^ (0) is given by fj43|) 
with 0" = /c (it is zero when (!32l) is imposed). Frequently (l43l) is written in a different form 
which is based on the following observation: if in the expression for the four- dimensional 
Ricci scalar R 

R = g-f^g^^'Ra.p. = gap,,, [g'^'g^" - g'^^g^") - \gapng,u,p 



X 



{g'^'^g^^g^" - Sg^^g^^g^" + 2^"''^^V'' + ^g^^g^pg^" - 4^"^^%'^^) , (44) 

we keep only the spatial indices and change the covariant component of g^"^ to e^"^ or, 
equivalently, impose the conditions fl32l) . we obtain the expression shown in square brackets 
of (H3ll . which is often called -R(3). 

Equation (H3i) gives contributions to the secondary constraints of zeroth order in the 
momenta p'^"^. There are obviously contributions to x^'^- Dirac's vector constraint, 7i^', does 
not have such contributions, so it is not directly related to the time development of the 
corresponding primary constraint p^^ (we will discuss this later). 

For (0), the equation (H3l) has to be compared to the corresponding expression of 
Dirac's (D39): 



Xl (0) = -B + {v^^7°°'''^?..,„i?™^^} ^ (45) 



where B (D38) 



B = \v^g''''"grs,ugab,v {i^-^^e™ + 2i?™'^^e^^} (46) 

is a part of full expression ([5]) where after passing to "e— form" only the terms cubic in e"^ are 
present. Terms quadratic and linear in e"^ are neglected, which results from the simplifying 
assumption because all non-cubic terms have either g^^ or the derivatives gokm- In his final 
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expressions ( H5|) Dirac keeps e*^™, not g''"^, which is consistent with his statement that this 
has to be true without the simphfying assumption which removes the difference between e'^" 
and g^"-. In addition, we keep g = det {gfj,u) in all equations. Dirac used = — det (g^u) 
and K"^ = — det (gkm) (or, probably, now more familiar notation '^g for det (g^^) and g (or 
(7(3)) for det (gkm) ) which are connected by (7°°J^ = or = a/— det ((7A:m) 7(7°° ■ 

By differentiating the second term of (H5|l . it is not difficult to derive the relation 

Dirac's scalar constraint Tii (Ti/, (0) = X/, (0)) is not the result of a direct calculation of 
{p^^,Hg}. This difference is not important for the proof of closure of the Dirac procedure 
and one can always consider linear combinations of constraints. For Castellani's procedure 
(or any other procedure) for finding gauge transformations we have to be careful with such 
redefinitions as we will demonstrate in the next Section. 

Until now, we have been concerned with the most complicated contributions to the sec- 
ondary constraints which are zeroth order in the momenta. Let us now consider the contri- 
butions to all orders. In the other two orders we obtain (using (l23ll ) 

X'" (2) = -^Hg (2) = i-L^ (gragst - Igrsgab) P'Y' , (48) 

X'"" (1) = ^^Hg (1) = -6: (^p:: - ie-V^^?..,. + e^^^.?™,.) . (49) 

Note, as (1) = there are no contributions linear in the momenta to the scalar 
constraint, but (2) 7^ 0, unless we impose ( l32i) . 

x''"" (0) was already calculated in fHSj) . For the full scalar constraint, x^^^ the relation 
(147|) is preserved in all orders 

The vector constraint has non-zero contributions in all orders of the momenta unless 
( l32l) is imposed. Before we continue to compare our direct calculations with those of Dirac, 
let us try to present the canonical Hamiltonian as a linear combination of the secondary 
constraints we calculated above. 
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We approach this problem by considering different orders in the momenta. The highest 
order is the second and the result is easily obtained from the first terms of fl23l) and fHHl) 



Hg (2) = (^graQsb - \grsga^ P'Y' = ^gOaX'" (2) (51) 

(using go^g^" = 5^ = 1). 

By considering fl27j) . which is equivalent to the second terms of fl23l) and fH9|) . we have in 
first order 

Hg (1) = -2guo (^p^/ - ^e"V^^/,,,, + e"V^^/™,,^ = 2go.x"'' (1) • (52) 

Hg (2) and Hg (1) are of the same form and we anticipate Hg (0) is also in this form. 
Unfortunately, this is not obvious and we have to perform some calculations to show it. To 
preserve the structure found in (ISTl [52l) . we will demonstrate that 



Hg (0) = -Lg (0) - Lx (0) = 2(7o.X°" (0) + (...),fc • (53) 

Note, that for Hg (1) given by fl27j) we also obtain fl52|) only up to a total spatial derivative. 
Using ([25]), ([38D, and (113]) we have 



-Lg (0) - Lx (0) - 2(7o.x"'^ (0) 



mnki „ / ~„ / „ukj j „viJ j 



-g^ gmn,i V i/i/Mi',* ['^ gOO ^ gOO 



(54) 



This equation demonstrates that the relations found for Hg (2) and Hg (1) are also valid 
for Hg (0) and the canonical Hamiltonian can be written in terms of x^'^as 

Hg = 2go^x"'- (55) 

Of course, this is correct up to total temporal (see ([H])) and spatial (see ([H]), ( 1271) . and ( IMIl ) 
derivatives. The modification of the initial Lagrangian ([8]) was proposed by Dirac while flM)) 
is obtained in the course of preserving relations found among contributions of higher order 
in the momenta to the constraints and the Hamiltonian. It would be very difficult to guess 
( I5^ without knowing the final result. Such an additional integration appearing in (|5^ . is 
very often performed at the Lagrangian level. For example, in the book by Gitman and 
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Tyutin [l7j, in addition to Dirac's ([8]) (which are B and first term of C* of Eq. (4.4.12) in 
171]), the integrations of fl5^ were performed at the Lagrangian level (the second and third 
terms of C*). The integrations of can be derived only in the course of the Hamiltonian 
procedure, but such integrations (if they are known) are also correct when applied to the 
Lagrangian because (going back to Dirac's derivation) it is clear that Lx (0) was constructed 
before the elimination of the velocities (i.e., at the Lagrangian level). 



How is this covariant form of Hq fl55l) (which is equivalent to what was found in 30l |) 
related to Dirac's expression for the Hamiltonian? Are they equivalent? The relationship 
between scalar constraints and Hl was found in (l50i) : we now consider the relation 
between the vector constraints. 

Let us inspect the form of our constraints calculated to different orders appearing in fl48l) . 
f l49|) . and fH3|) . There are simple relations between the contributions of different orders to 
and 

Ofc Ofc 

(2) = (2) , /°(1)=0, /nl)-V^^ X°nO) = (0) 

that allow one to write (to all orders) 



OA; 

X'" = + (56) 



with 



^Ofc ^ _pks _ ^kvprs Q^^^_^ _ g^^ ^^ (57) 

Solving fISB]) for tp^^ gives a combination of the constraints and which were originally 
calculated from the time development of the corresponding primary constraints. 

In terms of this combination of constraints ^^^^ and we obtain a different form of the 
canonical Hamiltonian 

Hg = 2^x'' + ^gok^'". (58) 
^uo 

This form of Hq is easy to compare with Dirac's, because his vector constraint is simply 
related to tp^'^ 

2^0fc _ gfc.-^^_ (gg) 
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For x"*^ we find 

Equation (l59ll . together with (H71) . demonstrates the equivalence of the two different forms 
of Hg given in ( l55i) and ( 158|) to Dirac's canonical Hamiltonian 

i^G = 2go„x''' = 2^X°° + 2^0^^°' = {g'') '^'^ + ^.oe^^K. (61) 

We would like to emphasize that Dirac's constraints are not a direct result of the time 
development of the primary constraints 0*^°" which produce ( (l50l) and ( l60l) ). The only 
place known to us where this is stated is in the book by Gitman and Tyutin (Eq. (4.4.19) 
of [l3]); but Dirac's particular combinations of constraints and the corresponding form of 
the Hamiltonian are usually used. 

The linear approximation of x^*^ gives exactly the constraints of linearized GR [37]. In 
the linearized case there is no difference between x"'^ and therefore linearized gravity 
can provide little "guidance" to full GR, in contrast to what was emphasized by ADM in 



4l|. Any such guidance has to be taken cautiously. 

To demonstrate closure of the Dirac procedure, any form of the canonical Hamiltonian 
( |6T|) is suitable as they are all equivalent; and any linear combination of constraints can be 
used for this purpose (e.g., = ^q'^'^^^^T^l and ip'^'^ = ^e^'^Hg)- When using Castellani's 
procedure to derive the gauge transformations generated by first-class constraints, we have 
to consider those secondary constraints that directly follow from the corresponding primary 
ones and the PBs of secondary constraints with the total Hamiltonian, not just with its 
canonical part (this is also discussed in the next Section). 

All of Dirac's secondary constraints have a zero PB with the primary constraints. In 
constraint dynamics this means that Lagrange multipliers cannot be found at this stage. As 
the PB of the secondary constraints with the canonical part of the Hamiltonian is zero or 
proportional to constraints, the procedure is closed. This is exactly the case here when we 
are taking into account the algebra^ of PBs among Dirac's combinations of the secondary 
constraints: 



^ This algebra is called "hypersurface deformation algebra" or "Dirac algebra" and can be found in many 
places, e.g. 0,0. 
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{Hl {x) , Hl (x')} = (x) Hs {x) {x - x') - e''' (x') K (a;') (x - x') , 

{x),Hl{x')}= Hl (x) 5 (x - x') , (62) 

{Tir (x) , Hs (x')} = Hs (x) (x - x) - Hr (x') 5 s(x') {x - x') . 

When dealing with the "covariant" secondary constraints x^'^ ^ the muhiphers are again 
not determined, but now we have 

{/^/n = ^^?'^V°. (63) 

The closure of Dirac's procedure is obviously preserved when using the covariant con- 
straints because x^"^ ^ind Dirac's constraints are simply related by fISUl) and fIBU]) . This can 
also be shown by direct calculation of {x^'^, Hq} without any reference to Dirac's combina- 
tions of constraints and their algebra. These calculations are long and to perform them we 
found it more convenient to work in the intermediate stages with x^^ and ip^^ . This allows 
us to sort out terms uniquely, and at the final stage we can express the result, using (l56l) . 



in terms of covariant constraints. The details of such calculations will be given in 421]. We 
arrive to the following PBs of x^'^ with the canonical part of the Hamiltonian 

{X ■j=lkmrbP QOad X + H ^9af3,kX ~ 7^9 900,bX (,04 j 

y 9 9 ^ 

and 



1 1 n'^^ 9 



{X ^ ^ i - -^J^-^ K^draP X - 9abP X ) ~ ~^~^J^^t^rbP 90ae X 



+ 9"'9oo,tx'' + 29op,t9V + ^9'' {9v,,r + 9rr,,, - 9r,,p) x'' - ^^?""^?oo,mX°°- (65) 

9 ^ 



Of course, we can present and (!65|) as one "covariant" equation 
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{X 7^ . ^tmrbP gOae X + '^m ; {^9raP X - 9abP X j+^oX,fc 

9 v~9 v~9 9 



- :,9'''9oo,bx'' + 9"'9oo,tx'' + 29op,t9"'x'' + ^9^' {9pq,r + 9rp,q - 9rq,p) (66) 
^ 9 

These equations, (jSH ESD or (jSSD, along with (jMD, provide proof of the closure of the 

Dirac procedure: higher order (tertiary) constraints do not appear and multipliers cannot 

be found because 



The Dirac Hamiltonian for GR, which is based on the modified Lagrangian of (fTOj) and 
the simplifying assumption fl32l) . is equivalent to the result of direct calculations given in fl55|) 
which are performed without any reference to surfaces of constant time. All the equivalent 
forms of the Hamiltonian of ( l6Ti) are only the consequence of an initial modification that 
does not affect the equations of motion and preserves the four- dimensional symmetry. It is 
natural to expect that Dirac's Hamiltonian formulation, which is obtained without any a 
priori assumptions and restrictions (e.g. surfaces of constant time), has to preserve another 
manifestation of four- dimensional symmetry: invariance under the diffeomorphism transfor- 
mation ([T]). Such a demonstration, given in the next Section, is an important consistency 
check of our results. All constraints are first-class and thus to find the generators of the 
gauge transformation, we have to consider "chains" of constraints. This means that one has 
to work, not with some combination of the constraints x^'^^ but with the exact results for 
and {x^'^jHt}- These results are complicated, especially (166|) . but their correct- 
ness can be verified if they lead to diffeomorphism invariance. A simple preliminary check 
of ( |66l) is that the linearized version of this equation gives 



{x'',H} = -p%^ {x'',H} = 0. 

hi 



This is equivalent to the results of [37| (note, that in the linearized case Xun = V'lirj.) it 
leads to a linearized version of diffeomorphism invariance. 

To summarize, the reality of Dirac's formulation, that is based on modifications of the 
initial Lagrangian, which do not affect the equations of motion, is as follows: notwithstanding 

27 



Dirac's references to space-like surfaces, all of his calculations were performed without use 
of any such surfaces. Consequently, Hawking's statement about the contradiction of the 
Hamiltonian formulation, based on splitting space-time into three spatial dimensions and 
one time dimension, is not applicable to Dirac's Hamiltonian formulation of GR, which does 
preserve the spirit of GR. Our own criticism of Dirac's formulation in {43! was not correct 
as we based it only on the 'interpretational' aspects of his work. This faux pas is also an 
illustration of how interpretations or some geometrical (or any other) reasonings can be 
dangerous if the "rule of procedure" referred by Lagrange is neglected. 

Dirac's simplifying assumption, gok = and ( l32l) . for constructing zeroth order in mo- 
menta contributions to the secondary constraints is not correct with respect to the individual 
parts given in (15^15^ : but remarkably when these parts are combined together in fH51) . they 
are equivalent to his final expression. His secondary constraints do not follow directly from 
the time development of the primary constraints but rather they are particular combinations 
of the true secondary constraints x^'^ ■ His secondary constraints cannot be directly used to 
find gauge transformations (Dirac did not consider himself this question). In the next Section 
we will show that the generator built from the true constraints gives the four- dimensional 
difTeomorphism ([1]), and we can say that the true constraints of the Dirac formulation and 
their algebra is "the algebra of four diffeomorphisms" 29 1. 



III. THE GAUGE GENERATOR AND TRANSFORMATION OF THE METRIC 
TENSOR 



The knowledge of the complete set of first-class constraints (primary, p^'^iQ, and sec- 
ondary, x^"" = x^"^ (2) + X^'^ (1) + X^'^ (0)) where contributions of different order in momenta 
are given by fHSl) . fH^ and fH5]) ). as well as the PBs between the primary and secondary 
constraints fl63|) . and the exact form of the closure fl66|) are sufficient to find the generators 
of the gauge transformations. This possibility is Dirac's old conjecture [3] which became a 



well developed algorithm and exists in a few variations 



11 



12| . We follow the work of 



id ] where the first application of such a method to Yang- Mills theory and ADM 



Castellani 
gravity^ was considered. 



10| the author referred to the Dirac formulation of GR but in fact considered the ADM formulation. 



In 

The non-equivalence of these two formulations will be discussed in the next Section. 
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Castellani's procedure is based on a derivation of the generator of gauge transformations 
which is defined by chains of first-class constraints. One starts with primary first-class 
constraint (s), i = 1,2,..., and construct the chain(s) ^•"■'G'*^^^ where ^j"'' is the nth order time 
derivative of the gauge parameter [n = 0, 1, ...). The maximum value of n corresponds 
to the length of the chain (e.g., n = 0, 1,2 for the system with tertiary constraints). The 
number of gauge parameters is equal to the number of first-class primary constraints. 
Note, that these chains are an unambiguous construction once the primary constraints are 
defined; the remaining members of the chain are uniquely determined. 

From this point, we specialize to the Dirac Hamiltonian formulation of GR with n = 0, 1 
and i = 0,1, {d— 1); if d = 4 there are four primary and four secondary constraints.^ The 
functions G^^^-j are calculated as follows 

G^i) (x) = p"'^ (x) , (67) 

Glo) i^) = + {x),Ht] + j aZ^ (x, y) {y) d^y (68) 

where the functions aZ^ (x, y) have to be chosen in such a way that the chain beginning with 
G^^j in ( 1671) ends on the primary constraint surface 

{GJq^, Ht} = primary. (69) 

The generator G {C,a) is given by 

G (^<t) = ^(tGJq^ + ^afiGJiy (70) 

There are some peculiarities that arise when applying this algorithm to GR that cannot be 
seen in simpler cases like Maxwell, Yang-Mills or linearized GR theories. 

Firstly, we comment on the use of different linear combinations of constraints. The term 
IpOo- ^ jj^^ ( |68ll is uniquely defined by the choice of primary constraints. After Dirac's 
modification ([8]) of the original Lagrangian, these remain just the momenta p^^ conjugate 
to the Qocr components of the metric tensor. Direct calculation of the PBs of these primary 
constraints with the Hamiltonian gives 

^ The following calculations, as well as the results of the previous Section, are valid in all dimensions, except 
d= 2. 
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Ok 1 -| 

r^Ofc /„\ 7T \ _ ^Ok _ ,nk I y ^ no _ ^„ksnj I /'„00\-V2 ^Ofc-1/ 

tP ia;j,/iG|-X - ~ 2 2 ^ 

and these expressions, X*"^? must be used when the gauge generators are derived. Of 

course, one can use Dirac's combinations, Tii and Tis, but only with appropriate coefficients 
or in appropriate combinations because of these inequahties: ^ ^l, {p^'',Hg} 

e^^Tis- We are not aware of any other situation where one must consider combinations 
of secondary constraints in ordinary field theories. Such a situation does not appear, for 
example, when the gauge generator for Yang-Mills is constructed lo|]. In this case, the ffist 
term of ( l68l) is just a secondary constraint which is the result of direct calculation of the 
PB of primary constraints with the Hamiltonian. In the Hamiltonian formulation of GR it 
is quite common (if not exclusive) to use of the Dirac combinations of constraints; but one 
has to be careful when gauge generators are constructed using T-Ll and Tis- 

Secondly, in both the Maxwell and Yang-Mills theories it is possible to choose the func- 



tions a (x, y) so that chains truly end with zero in flB^ 10|. This is not the case for GR, and 
chains end only on the surface of the primary constraints. The effect of such a difference 
will be seen in our calculation of the gauge generators and the associated transformations. 

Thirdly, the total Hamiltonian should be used in Castellani's procedure, not just its 
canonical part (see (!68l [69l) ). Again, in linearized GR, Yang-Mills and Maxwell theories 
this difference is irrelevant because in these theories the PBs of secondary constraints with 
primary ones are zero. This is not the case for full GR as can be seen from f l63p and similarly 



from equation (16) of [30 1. 



Finally, a purely technical comment. There is a change of sign in front of the ffist term of 
I 

(168|) relative to that used in \\S1\- This is the result of Dirac's convention for the fundamental 
brackets in ( !3TI) (it is the negative of the fundamental brackets used in lol|). 

To construct the generator (!70l) we have to find functions a'^ (x, y) using the condition 

{GI.^Ht] = j/'^ (x) + j a;{x,y)p''^ (y) d\ Ht 
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{x°'^(x),ifr} + j {a;{x,y),HT}p''^y)d'y + j a; {x,y) {p^^y) , Ht} d'y. (71) 

Part of the first term has aheady been calculated and (x) , Hq} is given by fl66|) . For 
the part involving the primary constraints, p^'^, we use (1631) which gives 

ix'" (x) , 9oo,op'' + 2gom,op''^ + He} = ^^7oo,o^7°'^X°° + 9omfl9'''^x'' 

^ T rJ^'^r, o'^^'-N/Or I ro- (r) rrflm Or ;ri«^-,/Om\ , rcr Ofc 

JtmrfeP fi'Oae X + r^ ~^ [^QraP X - QabP X j + ^0 



1 



The second term of flTTl) is irrelevant because it is automatically zero on the surface of 
primary constraints, as required by (!69l) . 

For the last term of ( ffTI) . taking into account the zero value of the PB among primary 
constraints, we find that 

{p''\Hr} = {p'\Ha}=x''. (73) 

This illustrates the advantage of using the 'covariant' constraints for deriving the gener- 
ator. With (I72H73I) . we can now read off the functions {x,y) from (I7T1) that compensates 



-< {x, y) = \9mfl {x) 9'" (x) 5; (x, y) + 9omfl9'''^5'^ 



Oc 2 11 



1 

- 1:9"'' 900 A + ^?°''^?00,i5; + 2^0p,i^^"5; + ^^7'^^ (^?p5,r + 9rp^, - 9r,,p) S',, (74) 

where 5^ (x, y) = 5^5 (x — y) and the arguments x and y are explicitly written only in the 
first term. Contracting (x, y) with the primary constraints and performing integration in 
the second term of fl68l) we obtain 
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^{0) — X - 2^00,09 p - gomflg p 



1 

+ :^^7'^'^?oo,;,P°° - g'^'gm^tp'' - 2gop,tg''V - ^g"' {gp,,r + grp,, - gr,,p) (75) 

g 

Equation (!75|) completes the calculation of the generator (!70|) . Now the transformation of 
fields can be found by calculating their PB with the generator^*^ 



S (field) = {field, G} . (76) 

For the time-time component of the metric tensor, qqo, we obtain (using Dirac's convention 
of (|3T|) and keeping only part of the generator flTOj) with terms proportional to p°°) 

Sgoo = {goo: G} = ~j^G = 
-Ku\g^Q<A-g p + gom,og p - -^g goo,bP ]+t,o,oP 



, _ _0m \ t , / l^_„^fcO , _ ^km\ p _ \j^h _ \ 



(77) 

Let us compare this result with diffeomorphism invariance ([1]) which can be written in 
an equivalent form, that is more convenient for comparison with our calculations: 

^(dijf)giiv = —^fj.,u — iv,^i + 5'"'^ {gtip,v + gup,n — gfj.i',13) ^a- (78) 

Taking /x = z/ = and explicitly separating space and time indices, we have 



Some authors defined transformations as 6 (field) = {G, field} which seems to be more natural. However, 
we keep the convention of Castellani that, of course, affects only an overall sign in the final result, which 
can always be incorporated into the gauge parameters (this is not a field dependent redefinition, as is used 
in 0). 
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S(diff)goo — —2.^0,0+ {g^^ 900,0 + '^g^''gok,o — g^^goo,k) ^o + {g^^goo,o + '^g^"^goin,o — g^'^goo,m) ik- 

(79) 

We see that ( j79l) is equivalent to (j77|) up to a numerical factor 2, 



2{5'oo,G'} = (5(rfi//)5(oo, (80) 

that can be incorporated into the gauge parameter by a rescaling ^„ 2^^. 
Similarly, for the space-time components, gok, we have 

Sgok = {gok, G} = 



5p 



Ok 



Om 



g"" \/-g 



-gg 



00 



1. 1. 



1 

2^ 



^"'^ 2 



1 1 



-g r 



Op 



+5' '^fi'OO.fc + "^gov.kg^" + Zm^""^ igpg,k + fl'fcp,g — fi'fcg,p) 



g 



(81) 



There is a difference between the transformation (|8T|) and the transformation of the 
time-time component (177|) as the momenta p"'' are present in (IHTl) . Using the definition of 
momenta f|T8|) and re-expressing e'^'" in terms of g^"^ by f|T5|) . we obtain 



+ [^°°^00,fc + (^Om,fc + 
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+ 2 b"^" (fi'OO.fc + gkOfi — gOk,o) + g"^" {gOn,k + gkn,0 — gOk,n)] C,m, (82) 

which again equals S(^diff)gfiu as given in flTHj) with fiu = Ok; which is true provided we again 
rescale by a factor of two. 

The last transformation to be checked is the transformation of the space-space compo- 
nents 



Sgkm = {gkm, G} = -^G ipP^) . (83) 



5 
5p 

The relevant part of the generator flTOj) which has an explicit dependence on jf^ is 



,Ocr 

G{f')=i.' 



g^^ V-g \f-g g 



(84) 

where parts of the secondary constraints (2), (2) and (1)) will also contribute to 
the final result. The variation of the last two terms in ( !84l) gives contributions proportional 
to the primary constraints (which equal zero on the surface of primary constraints). The 
only relevant parts of the generator for ^igkm are given by 

G = eoX°°(2)+a(x°'(2)+X°'(l)). (85) 

Performing variation of (!85|) with respect to p^"^ ^ using the expression for momenta given in 
f fTSj) . and reverting from e^"^ to g^"^ using f|T5l) . we obtain 



^gkm 2 (.^k,rn ~l~ ^m,fc) 



+ 2^*^° (fi'fcO.m + gmO,k — fi'fcm,o) ^0 + ^3^^ (fi'feO.m + fi'mCfc ~ gkmfi) 



+ 2^°^ (fl'fcp,m + S'mp.fc — gkm,p) ^0 + "fl'^'' (fi'fcg.m + fl'mg.fc " S'fcm.g) 'Cp- (86) 

It is not difficult to check that up to the same numerical factor 2 (as occurred in ( 180|) ) this 
is equivalent to 5(^diff)gp.u (CHI) with yuz/ = km. 

We see that transformations of the time-time and space-time components of the metric 
tensor are exactly equivalent to a diffeomorphism and the space-space components give a 
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diffeomorphism only on the surface of primary constraints. Such a deviation from ordinary 
field theories like Yang-Mills can be expected because the derivation of generators is per- 
formed (i.e., the functions aZ^ (x, y) are found) only on a surface of primary constraints. This 
is a consequence of the peculiarities of diffeomorphism transformations that will be discussed 
at the end of this Section. 

Returning to Dirac's statement [5] about abandoning four-dimensional symmetry in his 
approach; we can see that it is restrictive and only related to his initial modification of the 
Lagrangian ([8l). This abandoning of four- dimensional symmetry does not happen, neither 
in linearized [37| nor in full GR 30[. The Dirac Hamiltonian formulation of GR, as we 
demonstrated in this Section, allows one to derive the transformation of the metric tensor in 
covariant form and four-dimensional symmetry is preserved. The exact meaning of common 
statements, such as the one that is found in "unfortunately, the canonical treatment 
breaks the symmetry between space and time in general relativity", must also be clarified 
in light of our results. Of course, and this is a property of the Hamiltonian approach itself, 
the four-dimensional symmetry is not broken, it just is not manifest. For any generally 
covariant theory with first-class constraints we can only make a conclusion about abandoning 
such a symmetry if the gauge invariance that is derived from the first-class constraints 
cannot be presented in covariant form. This will be shown to happen in the case of the 
ADM formulation. If the symmetry presented in the original Lagrangian disappears in the 
Hamiltonian formulation, then it should be considered as a very strong indication that there 
is a mistake in the formulation; it is not a problem with the initial Lagrangian or with the 
Hamilton-Dirac method. From this point of view, if the "canonical treatment breaks the 
symmetry" , then such a treatment is not canonical. 

Let us return to the derivation of the gauge transformations. Our derivation of the 
transformation was based on an application of Castellani's method \ 
two variations of it: one of them is based on the extended Hamiltonian 



of. There are at least 



11 



ISj , where all first- 



class constraints are included, and the other, [12], is based only on the total Hamiltonian (i.e., 
only primary first-class constraints are included) as in Castellani's case. The e quiv alence of 
the algorith ms |ll | and [l^ was discussed in |l^ and a comparison of methods U] and |l^ 
was made in ll|] . Primary constraints play a special role in all of these methods. The need 
to include multipliers associated with the primary constraints was also emphasized in [ill . [isl 
and their importance in gauge transformations was demonstrated by some simple examples. 
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In 12], the multipliers are also important elements of this method (see below). Recently, the 



method of [l^ was applied to the ADM Hamiltonian in [l9|], where the transformation of 
the metric tensor was derived and was shown to differ in form from ([T]). For completeness, 
we apply this method to the Dirac Hamiltonian formulation of GR to demonstrate that 
the derivation of the diffeomorphism transformation in the Hamiltonian approach is not 
an artifact of a particular procedure for finding a gauge generator. At the same time this 
demonstration will illustrate the equivalence between the two different methods described 
in and 
problem in [26|]. 



as well as the equivalence both of them to the Lagrangian treatment of this 



The total Hamiltonian is the starting point of the method outlined in 12] 



Ht = Hc + X, 



(87) 



For a system with only irreducible secondary first-class constraints and no tertiary con- 
straints (so that {0^^, He} = x^) the generator of gauge transformations is simply 

G = ri,4>^ + i,x' (88) 

with two sets of parameters, and S,^ (twice the number of primary constraints), which are 
related by (see Eq. (17) of [12]): 



{s)fi 



(89) 



Here PF^, V^, C^, and B'P' are the structure functions of the involutive algebra (see Eqs. 
(2) and (3) of [jj) 



(90) 



{i^e,x^} = v^(;;),0'^ + V(:),x^ 



(91) 



(92) 



(93) 
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The indices (p) and (s) indicate structure functions associated with the primary and 
secondary constraints respectively. Note, that fl89l) involves only functions with the 
subscript (s), i.e. structure functions related to the primary constraints are not present, so 
that this equation is valid on the primary constraint surface. This is similar to what happens 
in Castellani's procedure. 

To find the generator (!88|) . one has to solve (l89l) for rj^; and as in Castellani's method the 
number of independent parameters becomes equal to the number of primary constraints. For 
the Dirac Hamiltonian of GR, which is obtained after modification of the gamma-gamma 
part (with no effect on the equations of motion, canonical variables ^f^jy, and conjugate 
momenta p'^'^) we have the simple primary constraints 

/'^ = 

and the secondary constraints (for explicit expressions see (HHl) . (H9|l and (j43|l ) 

{0ri/,} = /^ (94) 

This allows us to write the Hamiltonian in a compact and symmetric form 

He = 2go,x'"'- 

The possibility of solving (189|) for i]^, which is an ordinary algebraic equation, depends on 
the structure functions Wy^ and Cfl . In the case of the Dirac Hamiltonian formulation 
of GR (see (EID) they are 

n), = ^;^ ^S. = (95) 

and we can solve (|89|) for r]^: 

V, = - (^(:), + A,5(7),) . (96) 

The structure functions V^^-j^ for GR are complicated (see ( 166|) ) and for -^(^^ we have (see 

m) 
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and (p6|) becomes 



V, = ^,,0 - ^uVC), - ^ugo.filg'-'Sl. (97) 

Note that the structure function B'^^)/! ^'^^^ equal zero and the Lagrange multiphers, 
which are the velocities goa,o that cannot be expressed in term of enter (1971) explicitly 
(■^0 = fi'oo.o and Xk = 2gQk Q)}^ After the substitution of fl97|) into flHHl) we obtain a one- 
parameter (the number of components equals to the number of primary constraints) 
generator 



G (6) = -6^?oo,o^^7^V° - ^u9om,o9''^'<P' + U^" - ^uVi;,)^ct>^ + e.X^ (98) 



This expression has to be compared to the generator found using Castellani's procedure 

G = ^i/,oG'('i) + ^lyG'^Q^ . (99) 

The third term of (198!) is exactly the same as the first term of (1991) : the first, second, and last 
terms of fl^Sl) are also the same as in Castellani's approach (see the first line of fl75]) ). The 
fourth term of ( |98l) . to be compared with our result obtained using Castellani's approach, is 
—(.uV^g-^^4>^. The structure function, V'^g)^, originates from the calculation of the PB of ( l9Til . 
which in the case of a field theory is 



{H,,x''{x)} = J V^:)^{x,y)x'{y)d'y. (100) 

The direct calculation of {He, x'^} gives terms proportional to x'^ and its derivatives; for its 
explicit form see (1661) 

{He,xl = K;x' + M;'x':k- (101) 

This equation is usually presented in the following form in order to find the structure func- 
tions 

{He, x^} = I (^K; (x) 6{x-y)- M;:' (x) ^5 (x - y)^ x' iv) d'y. (102) 



Expressions for the multipliers came from the Legendre transformation and from the fact that the modified 
Lagrangian is independent of ffo/^.o, so that H = gut^fiP"^ - i = ffoo.op"" + 2^0/^,0^"^ + •■• 
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This is the standard form for an intermediate result in such calculations (e.g. see |10|). 
To find the generator fl98l) . it is not necessary to rewrite (110 II) in the form of (I102p : and 
the direct substitution ~^ iiito (110 II) gives the corresponding part, —^uy(^s)ti^^y 
generator fl^S]) . This is equivalent to the second and third lines of fITS]) . 



If the novel method of [1^ is applied to the Dirac Hamiltonian of GR and not to the ADM 
Hamiltonian, as was done in 19(], then the generator (!98|) is equivalent to the one obtained 
using the old procedure of [lOj]; and so they generate the same gauge transformation. We 
have found that the results for the gauge transformation are equivalent, whatever method 
is used to find it, when applied to the same Hamiltonian. 

The peculiarities of Castellani's procedure, when it is applied to the Hamiltonian of GR, 
were discussed at the beginning of this Section and illustrated by a derivation of ([T]) using the 



methods [10| and [12]. They are originated from the algebra of constraints either in Dirac's 



formulation, given in this Section, or in the 'covariant' formulation of jSOj. This algebra is 
different from that of ordinary gauge theories. It reflects the peculiarities of diffeomorphism 
invariance if compared to the gauge invariance of ordinary gauge theories. We now briefly 
consider this topic. 

In the Introduction, we restrict our discussion to a particular meaning of diffeomorphism 
given in ([1]) that is generally accepted in literature on GR and which is similar to the usual 
gauge transformations. It is in exactly this sense that diffeomorphism invariance can be 
derived from the Hamiltonian approach. Now we would like to describe this transformation 
without recourse to the Hamiltonian formulation as it is usually presented in textbooks (e.g. 



see 201]); and we wish to reveal how a difference between these two views of diffeomorphism 
invariance manifests itself. This exercise will also demonstrate the connection between ([1]) 
and general coordinate transformations. 

The principle of general covariance, the cornerstone of GR, puts severe restrictions on 
the possible forms of the Lagrangian. The simplest is the EH Lagrangian.^^ The EH action 
and the Einstein equations are invariant under a general coordinate transformation 

^'f^ = (2;'') (103) 
and the corresponding transformation of the metric tensor 



Of course, it is not unique and there are many posibilities: such as Lovelock gravity 4J] or / (i?) 
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For infinitesimal transformations 

^ x'^ = x^' + e {x) (105) 

(11041) can be written as 

g"" {x') = g'" (x) + er.^?'^" {x) + er^^"'^ (a:) + O {e) ■ (106) 



Note tliat tlie components form a true vector j38l. I45l|. in contrast to tlie parameters e 
and which appear in the ADM formulation of j^J] (see ([2]) and its derivation in the next 
Section). 

If we consider (x) as being a small parameter and restrict our interest to the first-order 
contributions in then the exact invariance of the EH action and the Einstein equations of 
motion is lost as only the inclusion of all terms in the expansion of (I106P will preserve it. In 
addition, if we want to present (11061) in a form similar to a gauge transformation, in which 
the invariance with respect to replacement of variables is in the same coordinate frame of 



reference 461], we should write both sides of (11061) in the same coordinate system. This can 



be done by an additional approximation, using the Taylor expansion of g' {x'): 



g'i'^ {x') = g'"" (x" + C (x)) = g'^^ (x) + gf^^C + O {^^) (107) 
where in the second term (which is already linear in ^") we used ^^g^/^a ^ = g^a + ^ (0- 

x'=x ' 

Combining (I106p and (I107p and keeping only terms linear in we obtain 



Sg^-" = g'^^ (x) - g^^ (x) = C^?"" {x) + i^g^^ (x) - g^^C (108) 
which is equivalent to 

By repeating similar calculations, or by using b {gnag""^) = 0, one can find transformations 
for the covariant components of a metric tensor 
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Sg,u = -^,;u - = -^M,- - ^u,, + 2Vl^U (110) 

This equation is just ([T]) (or its equivalent form ( |78|) that was more convenient for comparison 
with our calculations). With being a true vector, the transformations (11091) and (1 11 01) are 
generally covariant (these are covariant derivatives of a true vector), so they are independent 
of the choice of coordinate system; these are the transformations we derived from the Dirac 
Hamiltonian of GR and in 30 1. 



Similarly, we can obtain the transformation of the Christoffel symbols. Using the relation 
between V^^ and g^j_y in (fT3l) and the transformation 5gij_y in (1 11 01) . we obtain 



Note that the presence of second-order derivatives of the parameters (^' 



'1111 



im- 



mediately allows one to come to a general conclusion about the constraint structure of the 
Hamiltonian formulation of GR in the first-order form, the Einstein affine-metric formulation 
of a. where .-^ and PJ^ a.e Seated a. independent fields The presence of the seeond- 
order derivatives of the parameters in the transformation SV^^, implies that the generators 
must have the same order of derivatives, i.e. the tertiary constraints must appear in such a 



formulation. Of course, direct calculations confirmed this simple observation 
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51 



52 



53|. 



In the first discussion of the Hamiltonian formulation of the first-order form of GR which 



was presented in 4l|, the tertiary constraints did not appear because some first-class con- 
straints were solved before closure of Dirac's procedure was reached; this procedure is not a 
consistent implementation of Dirac's procedure for treating constrained systems. This fact 
was clearly demonstrated in the Appendix of 5J] and was discussed in 43|. The loss of 
tertiary constraints is also due to a misleading analo gy b etween the first-order formulations 



of Electrodynamics and linearized GR appearing in [41|, |55|. In the first-order formulation 
of Electrodynamics, where the field strength is treated as an independent variable, there is 
no increase in the order of the derivatives of the gauge parameters in the generator of gauge 
transformations because this is relates to the fact that the variation of the field strength is 
zero. In contrast, the variation of F^^ under diffeomorphism transformations, (lllll) . is not 



^•^ This formulation was originally introduced by Einstein js^ (for English translation see [l^]), but mistakenly 



attributed to Palatini 



48 1 (see also Palatini's original paper {49] and its English translation 50]). 
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zero and cannot even be written in covariant form as Sgf^i, was in flllOp . This characteristic is 



consistent with FJ^^, not being a true tensor 



20 



56| ■ Also there is no hnear combination of 



the first-order derivatives of the metric tensor that is exactly invariant under diffeomorphism 
transformations, as well as under general coordinate transformations.^^ This invariance is 
also related to the fact that a generally covariant action for GR cannot be built from terms 
only quadratic in the first-order derivatives of the metric tensor; and the simplest generally 
covariant, EH Lagrangian, is proportional to a Ricci scalar and this involves second-order 
derivatives of the metric tensor. This is why the affine-metric formulation of GR leads 
unavoidably to tertiary constraints and consequently increases the length of the chain of 
constraints and the order of the derivatives in the gauge generator. 

Let us compare fillip with the gauge invariance of Yang- Mills theory. We have the field 
strength F^^, whose variation does not vanish, in contrast to Electrodynamics, 

with /"''^ a totally antisymmetric structure constant and 6^ a gauge parameter. We do 
not have exact invariance for F^^, but the gauge parameters enter only linearly. Thus 
in the first-order formulation of Yang-Mills theory, if we consider the field strength as an 
independent variable, there is no increase in the length of the chains of constraints needed 
to accommodate these transformations. In GR it is not possible to build any combination 
of first-order derivatives which is exactly invariant under diffeomorphism; and it is also 
impossible to find a combination whose variation is proportional to the gauge parameter or 
its first-order derivatives. 

The transformation ([1]) (or other equivalent forms) is written in the same coordinate 
system and, because this combination is a true tensor, it is independent of the coordinate 
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58|; but 



transformations. In this sense it is "analogous to the gauge transformation" 
the absence of combinations of derivatives that do not lead to an increase of the order of 
the derivatives of gauge parameters in the generator is a distinct property of GR. 



The r" behaves hke a tensor only with respect to hnear coordinate trasformations 20|. Probably, this 



was the origin of the analogy between Electrodynamics and linearized GR and of the conjecture that this 



analogy should be extended to full GR [41 1 . 

Actually, such a combination exists which is a true tensor: this is the covariant derivative of the metric 
tensor, g^v;'y, but it identically equals zero. 
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In addition, the Lagrangian of GR is not exactly invariant under a diffeomorphism trans- 
formation, in contrast to Maxwell and Yang-Mills theories. From (11081) and fillip and by 
using 



D -pa "pa I "pa "p/3 -pa p/3 

we can find the transformations of and R = g"'^Rap'- 

5R^, = -eR^.u,p-^':,Rup-^'.R,p, SR=-CPR,,. (112) 

From the above relations it is not difficult to obtain the transformation of the EH Lagrangian 
under a diffeomorphism 



6 i^R) = i-ev^R) 



;ii3) 



This lack of exact invariance is distinct from what occurs in the Maxwell and Yang- 
Mills theories; but other models exist with Lagrangians which are also invariant up to a 
total divergence, e.g. Topologically Massive Electrodynamics (TME) of [59|. (See 
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621 ] for a discussion of its first-order formulation.) Despite differences in the invariance 



property of Lagrangians, which can be either exact as in ordinary Electrodynamics or up to 
a total divergence as in TME 59|] , the equations of motion are exactly invariant under gauge 
transformations. In GR the transformation of the equations of motion is proportional to the 



equations themselves 
field equations 



46| . Using (11121) we can easily find transformations of the Einstein 



SG,, = -eC^^.p - ^'pG^p - ^'uG,p, (114) 

where = R^u — \gp.vR is the Einstein tensor. So, in GR under diffeomorphism transfor- 
mations ([T]), the equations of motion are invariant only 'on-shell' which does not contradict 
the principle of gauge invariance: a solution to the equations of motion maps into a solution. 
This last result, the 'on-shell' invariance of the equations of motion might cause some con- 
fusion because ([1]) was obtained from infinitesimal coordinate transformations (llOSp . which 
are a particular case of general coordinate transformations (llOSp . The EH action and the 
Einstein equations are exactly invariant under (11031) . therefore, under (11051) (in fact, the 
Einstein equations were originally postulated so as to satisfy (11031) ). After writing (llOSp in 
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the same coordinate system and using the approximations (11061) and fll07p . exact invariance 
is lost. Such a mapping, (11141) . by itself, is not peculiar to GR because a similar property is 
present in the Yang-Mills theory where the transformation of the equations of motion is 



6iD,Fn=fabcetD,Fr. 

As in GR, the Yang-Mills field equation is only invariant 'on-shell'; but this transformation 
is proportional to the gauge parameter, whereas the transformation of SGf^^, in (11141) also 
contains derivatives of the gauge parameter. The peculiar algebra of constraints in GR is 
related to this increase in the order of the derivatives of the gauge parameters in the trans- 
formations of the equations of motion, as well as the impossibility of finding a combination 
of derivatives of the metric tensor which is either exactly invariant or whose variation does 
not require derivatives of the gauge parameters. 

On one hand, the diffeomorphism transformation ([1]) is related to coordinate transfor- 
mations and, as a consequence, to general covariance. On the other hand, the Hamiltonian 
formulation of GR is performed using the same "rule of procedure" as in ordinary gauge 
theories, allowing one to derive the same transformation ([T]) . The resulting invariance makes 
this transformations distinct from ordinary gauge theories because of the presence of deriva- 
tives of the gauge parameters. Such distinct transformations should manifest themselves in 
all steps of the procedure and this is exactly what we have found. For example, the non-zero 
PB among primary and secondary constraints in ( l63|) and the 'on-shell' of primary con- 
straints result of (18^ are properties that are not observed in ordinary gauge theories. These 
peculiarities are present in Dirac's formulation as considered in this article, as well as in 30|; 
and both of these formulations allow one to derive the diffeomorphism invariance ([1]). At 
least two of the above mentioned properties, which are related to (1631) and (l84l) . are absent 
in the ADM formulation and the transformations derived from the ADM Hamiltonian are 
different from a diffeomorphism. A comparison between the Dirac and ADM formulations 
will be made in the next Section. 
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IV. THE DIRAC HAMILTONIAN OF GR VERSUS THE ADM HAMILTONIAN 
OF GEOMETRODYNAMICS 



In the previous Section we demonstrated that, following the "rule of procedure" as ap- 
plied to Dirac's Hamiltonian of GR flMl) (with the covariant metric g^i, and corresponding 
conjugate momentum p^'^ as the fundamental canonical variables) the diffeomorphism in- 



variance ([T]) can be derived using Castellani's procedure 10 1 or the method of ll^- The 
same result was recently obtained in the Hamiltonian formulation of GR [30| without using 
Dirac's modifications (IHl) and in the Lagrangian approach of by Samanta {26]. This 
result is exactly what one would expect for the invariance of GR, as well as the equivalence 
of results in the Hamiltonian and Lagrangian approaches. 

In addition, because Dirac's non-covariant modifications of the Lagrangian do not change 
the equations of motion, four- dimensional symmetry is preserved and it is reflected in the 
covariant form of the transformations ([1]). We have also demonstrated that Dirac's refer- 
ences to space- like surfaces are not part of his actual calculations. As a result, Hawking's 



statement 27|] that "the split into three spatial dimensions and one time dimension seems 
to be contrary to the whole spirit of relativity" is not related to Dirac's formulation, where 
only manifest invariance (but not the invariance itself) is broken by explicitly considering 
different components of the metric tensor. We have seen that working with the original 
Einstein variable, the metric tensor, we have the Hamiltonian formulation of GR that is 
consistent with diffeomorphism symmetry, and the spirit of GR is "alive" . 

There exists a more popular Hamiltonian which is based on a different set of variables: 
the lapse and shift functions and the space-space components of the metric tensor. It is 



frequently, but mistakenly called the Dirac Hamiltonian (e.g. 



idp and even its variables are 



called "Dirac's lapse and shift" [22]. This formulation (with the lapse and shift functions) 



is due to Arnowitt, Deser and Misner (see jo] and references therein). The name "Dirac- 
ADM" is also not correct; and despite the apparent similarities between the Dirac and ADM 
Hamiltonians, they are different (see below). The appropriate and best known names for 
the ADM formulation are "ADM gravity" and geometrodynamics, as opposed to Einstein 
gravity. The gauge transformation derived from the ADM Hamiltonian by the methods of 



16 



The names "lapse" and "shift" were introduced neither by ADM nor by Dirac and appeared only later, 



for example, in [l4l. |32||. 
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id . Il2j, using ADM variables and constraints, is not diffeomorphism invariance ([T]). This 



fact was recently demonstrated very clearly using the method of [1^ in [19| and was discussed 
in our Introduction. The field- dependent redefinition of the parameters ([2]) for the ADM 
Hamiltonian differs from the exact result that was obtained in Section III from the Dirac 
formulation, where no field-dependent redefinition was used, and in |2Q] and [s^ where the 
question of equivalence with diffeomorphism does not arise. 



We feel that it is insufficient to say that the formulations of 14, 1301 and that of Dirac 
[5|], both of which use the metric tensor as the canonical variable, correctly describe the 
Hamiltonian of GR and one is obliged to work with the original Einstein variables. One 
should not try to recast GR into a description of the motion of surfaces; more precisely, 
one should not change variables to make such an interpretation plausible. It is necessary 
to understand why two such closely related Hamiltonians, these of Dirac and ADM, which 
are mistakenly said to be equivalent, produce different results. This is a general question 
related to the Hamiltonian formulation of singular Lagrangians. The understanding of the 
peculiarities of the Hamiltonian method for constrained systems can protect from repeating 
some mistakes that have been made when considering formulations of such theories. 

The analysis of the ADM formulation is also interesting from another point of view, 
as it provides a very instructive example of what might be called an "interpretational" 
approach to physics. In the original work of ADM ( e.g., [^), an interpretation of the 
variables they introduced was proposed; and in later work, this interpretation became the 
cornerstone of that theory. Attempts were made to "derive" results starting from that 
interpretation, i.e. by elevating the interpretation to the level of first principles 321]. Such 
an approach exhibits lack of rigour and relies on some geometrical reasoning. The essence of 



this approach is perfectly expressed in the following quotation from [63| : "I capture as much 
of the classical theory as I can by pictorial visualization" ^ and "The reader is encouraged to 
follow the broad outlines and not worry about technical details" . The "advantage" of such 
an approach is that it cannot be disproved; yet it prevents one from obtaining any reliable 
prediction or result. For example, if we accept Dirac's references to space-like surfaces as 
a part of his formalism, then Hawking's statement that introduction of a family of space- 



This becomes a well-known pedagogical approach in teaching conceptional (not mathematical) physics for 
non-science students. 
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like surfaces "seems to be contrary to the whole spirit of relativity" 27| forces us to reject 
this formulation as an inappropriate formulation of GR. And yet, as we have demonstrated, 
Dirac's formulation leads to a direct restoration of diffeomorphism invariance and, because of 
this, it is consistent with the spirit of GR and is the correct Hamiltonian formulation of GR. 
This demonstration shows how a purely interpretational consideration can lead to a wrong 
conclusion and that the interpretational approach without having to "worry about technical 
details" is meaningless. Dirac's derivation follows the "rule of procedure" and allows us 
to check any interpretation by explicit calculations. If something is constructed only using 
pure interpretation, then the final result cannot be checked by calculations and can only be 
analyzed by comparing its consistency with general principles. An "interpretation" cannot 
serve as a ground to disprove a result and in fact could be wrong, as in the case of Dirac's 
formulation in which he makes references to space-like surfaces that were not used in his 
calculations. A general understanding of the limitations of the interpretational approach 
probably provides an explanation of why Hawking's words 23], spoken almost thirty years 
ago on the occasion of the centenary of Einstein's birth, were not enough to cause people to 
immediately abandon the ADM formulation. 

Yet more disturbing, is that this "interpretational" language has completely prevailed in 
the Hamiltonian formulation of GR. As an example, consider the correct work of Samanta 



261] in which he used the Lagrangian formulation, and where there are no surfaces of constant 



time, space-like, slicing, etc., then it is abruptly altered when he refers to the Hamiltonian 
formulation of the same theory and states that "slicing is essential for Hamiltonian formu- 
lation". This assertion is obviously wrong as slicing is not essential and the Hamiltonian 
formulations of GR obtained without any reference to slicing gives a consistent result, as we 



30]. 



have demonstrated in the previous Section and in 

The general trend in Physics and the main goal of many physicists is the unification of 
theories and methods on all possible levels; but even now, when we are a few years away 
from the 100th jubilee of the discovery of GR 64] there remains an inconsistency when 
discussing of different formulations (Lagrange and Hamilton-Dirac) of the same theory, Ein- 
stein's GR! This leads to erroneous observations, such as "It is worth noting that generalized 
Hamiltonian dynamics is not completely equivalent to Lagrangian formulation of the origi- 
nal theory. In the Hamiltonian formalism the constraints generate transformations of phase 
space variables, however, the group of these transformations does not have to be equivalent 
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to the group of gauge transformations of Lagrangian theory" [65 1. 

We consider it important to understand and find exphcitly where and why the ADM 
Hamiltonian formulation contradicts the spirit of GR, and why it cannot be associated with 
Einstein's theory (i.e. it is not the Hamihonian formulation of GR, but rather the Hamilto- 
nian formulation of distinct theory: "geometrodynamics"). One can talk about abandoning 
he spirit of GR if one is discussing a different theory built on different principles (such as 



661]) but the ADM formulation has been given the appearance of having a formal basis on 
GR (in view of their articles and the presentation in many textbooks); indeed the summary 
of ADM's work in has the title "Dynamics of General Relativity". The mathematical 
manifestation of the spirit of Einstein's GR is the general covariance of Einstein's equations 
of motion for the metric tensor. Einstein's GR is a field theory and the methods used in or- 
dinary field theories, those of Lagrange and Hamilton, if applied correctly, must not destroy 
its spirit. This is exactly what we want to investigate: where was "a regular and uniform 
rule of procedure" broken in the ADM approach. In Hamiltonian language, we want to see 
where the canonical procedure was destroyed by passing from the Dirac Hamiltonian to the 
ADM Hamiltonian and why the ADM formulation with their variables is not equivalent to 
GR or, in other words, is not a canonical formulation of GR. 

It does not seem possible to start from the Lagrangian of GR, where surfaces are not 
present, and then after introducing new variables have such surfaces. We will follow a path 
distinct from the interpretational approach and pay attention to technical details by using 
the "uniform rule of procedure" in analyzing the ADM formulation. In the previous Section 
we demonstrated that, when we are using the rule of procedure, surfaces do not appear in 



either Dirac's calculations or in [30|. We will not repeat the calculations of the previous 



30l |. but instead we will analyze how the ADM formulation is related to that 



Section or of 
of Dirac. 

Let us compare the two Hamiltonians of Dirac and ADM. Castellani himself considered 



the GR and Yang- Mills theories as an illustrative examples o 



transformations. Referring to Dirac's book [3], Castellani 1^ started with the statement 
"from the Hilbert action one derives the Hamiltonian" 



his algorithm for finding gauge 



Such a Hamiltonian can be "derived" without recourse to the EH Lagrangian. We refer the reader 
interested in "visuahzation" or in the geometrical reasoning behind a derivation or geometrical meaning 



of this equation to the numerous figures in S^]. We disregard such approaches as inadequate for any 
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HAUM ^ N^AUM ^ ^^^AVM ^ ^^^^ + N^^Ui, (115) 

where 11 and Ilj are momenta conjugate to and N^, respectively. 



This Hamiltonian never appears in Dirac's 



of the Hamiltonian of GR is in the article js] that we discussed in previous Sections. It is 
different from flllSp and given by fl6T]) which is the canonical part of the Hamiltonian, and his 
primary constraints are the momenta p^^ conjugate to go^ IQ. Equation (11151) is, in fact, the 
ADM result and we have indicated so by using the superscript ^ADM\ In order to compare 
the results of the Dirac and ADM formulations we use a slightly different convention for the 

n 

Dirac Hamiltonian which appeared in a subsequent article of Dirac (see Eq. (7) of [67] for 
the canonical part of Hamiltonian) 

Hd = i-g'') '^'^ Ul + Qroe^'Hs + <7oo,oP°° + 2^7ofe,oP°'- (116) 



Dook or in any of his articles. Dirac's derivation 



In this convention, qqq is negative [6 
Let us compare the secondary co 
morphism" constraint, is given in many sources (e.g. see Eq. (3.14b) of [6.]) as 



Let us compare the secondary constraints of the two formulations. 7i^^^ , the "diffeo- 



where H'^-' is a momentum conjugate to the spatial metric g^j- The symbol "|" seems to 
indicate the covariant derivative with respect to g^j 6|; but the definition of the particular 
covariant derivative used in ADM is non-standard and is not easily found (see Eq. (5.1) of 



68( 1 and Eq. (2.3b) of 69| which are consequently papers six and twelve in their series) 



yrkj _ T-rfcj I ytlm-pk 

\j — ,j Im 

which gives 



This is exactly the Dirac constraint (see Eq. (D41) of [Sj]) or our (1291) ). We note that the 
ADM definition of a covariant derivative of a contravariant second rank tensor mimics a 

proofs. 
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covariant derivative of a contravariant vector (the first rank tensor) or a covariant derivative 



of the tensor densit y [38 
second rank tensor uQ 



38 



f we use the standard definition of a covariant derivative of the 



56| we will obtain a different result. This result can only be 



presented as a standard covariant derivative (but with respect to the spatial metric g^k only) 
if we write it as 



n 



Dirac _ '^ADM 



(see Eq. (37) of [70] or Eq. (E.2.34) of (l5|). 



Vdet 



9k 



The scalar, "Hamiltonian" , constraint TC^^^'^ is given by Eq. (3.14b) of Gj 



n 



ADM 



-V9 



Vg^R + -^[ gik9jm - 7^9ij9km ] wm''"'. 



;ii7) 



Using (50) (employing the convention of 67| where g^^ is negative), from (48) we obtain 
the second term of flll7p . Also from (43), and taking into account (44), the first term of 
(11171) follows. Thus, this constraint is also equivalent to Dirac's Ti^. 

Dirac's combinations Hl and Hs of the true secondary constraints (given in (l50 l [60l) ) 
are exactly the same as the ADM secondary^^ constraints TCf^^^ and Tif^^'^ 



H 



ADM 



ADM 



Hi 



(118) 



The only difference between the first two terms of flll5p and flll6p is that the field- dependent 
coefficients in front of Dirac combinations of constraints are called new variables by ADM. 
These are the lapse and shift functions 



00\-l/2 



;ii9) 



9 



Oi 



(120) 



Actually, in the ADM formulation they appear as primary 0, [t^ • We will return to this later. 
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In addition, Dirac's e-'* (fT5l) (the reciprocal of gji) is called the "three-dimensional" metric 
g^^ j^. (In fact, g^^ is the space-space component of the full four-dimensional metric in four- 
dimensional space-time.) To distinguish e-'* from the space-space components of Einstein's 
four-dimensional contravariant metric tensor, the latter is defined to be '^g^^ Q]. Dirac's 
notation is more transparent as it arises in his derivation of the Hamiltonian that was 
analyzed in Section 11.^'' 

The confused notion (in the literature) that these two formulations are equivalent, is 
understandable, especially in light of relation (11181) . Another reason is that in many presen- 
tations of the ADM Hamiltonian (11151) . such as in (29I, the primary constraints are ignored by 
imposing the idea that the lapse and shift variables are merely the Lagrange multipliers for 
the constraints 7i_L and Tij and that they can be treated as 'primary' rather than 'secondary' 
constraints. Firstly, in such an approach even derivation of the gauge transformations of all 
components of the metric tensor becomes impossible as we are no longer dealing with the 
full phase space of the Hamiltonian (e.g., see the remark 011 p. 3288 of [2^). The methods 
of derivation of Castellani [l^ and of 12 1 (which is used in 19 1) employ the complete phase 



space and so all fields and their conjugate momenta must be included. Secondly, dropping 
primary constraints contradicts the methods of constraint dynamics: primary constraints 
are first-class and must not be solved as this destroys the gauge invariance present in the 
Lagrangian. Only second-class constraints can be solved and then some of the variables 
can be eliminated provided PBs are replaced by Dirac brackets. All this means that if we 
derive generators of gauge transformations, following any procedure using Tif^^'^ and Tif^^ 
as the first-class primary constraints, we will have generators independent of the momenta 
conjugate to the "multipliers", so that the gauge transformation of, for example A^, would 
be zero 

SN = {N,G} = 0. (121) 

This result means that 6g^^ = 0, which is not related to a diffeomorphism transformation 
at all. Without the primary constraints the time derivatives of the lapse function, for 
example, would vanish according to the Hamiltonian formulation. 



The ADM renaming was probably introduced to underline the geometrical interpretation of their variables. 
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{N,H} = 0. 

Thus, is constant in time, yet if we use the total Hamiltonian (11161) . we obtain 



{N,Ht} = N,o. 

Finally, in the Hamiltonian approach, the primary constraints come from varying the 
Lagrangian with respect to velocities, and if we follow this rule, then Ti^-^*^ and Ti.f^'^^ are 
not primary constraints. We will consequently work with the total Hamiltonian. 

Let us continue to compare these two formulations. The relation between their primary 
constraints is as yet unclear and we shall return to this later. Form (11191) and (11201) the 



metric and its inverse are (e.g., see 



3,Q 



9 



1 



Af2 



Af2 y Ar2 



(122) 



9ijN^ 9ij 

The fundamental PB of the canonical variables, the components of the covariant metric 

n 

tensor and their corresponding conjugate momenta, for Dirac's Hamiltonian are |5| 



(x) , g,^ {x') } = ^ (5^5,^ + 5X) 5, {x - x' 
whereas for the ADM approach the fundamental PBs are (e.g., see 



10 



19 



nl) 



{g,, (x) , H^' {x')] = \ {5^5\ + 6, {x - x') = ^5, {x - x') 



(123) 



(124) 



[N' (x),H, (x')} = 5]5^{x-x') 



(125) 



{N{x),Ii{x')} = 5^{x-x'). (126) 

Other PBs presumably equal zero (i.e. (x) , H'^^ (a;')} = 0, etc.) if these variables are to 
be canonical. 

We now investigate why the Dirac and ADM approaches to the canonical treatment of 
GR lead to diffeomorphism transformations in the former case and to transformations that 
correspond to a diffeomorphism only after a non-covariant field-dependent redefinition of 
gauge parameters in the latter case. 
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In Section III we have demonstrated that the gauge transformations that follow from 



Dirac's Hamiltonian can be derived using both the methods of and 12|. The method 



12 



n 

was applied to ADM gravity in [19!']. The result of [I9l] is not new, but it is probably the first 
complete consideration of how one can derive the gauge transformations from the constraint 
structure of the ADM formulation. The application of Castellani's method to the ADM 
Hamiltonian given in Appendix of [lo| is opaque and incomplete as the relation between the 
diffeomorphism and the ADM parameters is not explicitly given and only the transformations 
of go^ are found. The calculations themselves were performed in an unnatural way - in order 
to find the transformations of the metric tensor, the ADM variables were expressed in terms 
of the metric in the generator. For completeness, we shall use Castellani's approach to find 



the gauge generator with the ADM Hamiltonian and compare it with the result of p^. In 
addition, we will show some of the peculiarities in this calculation which are related to the 
somewhat confusing notation used by ADM. 

According to Castellani's procedure, the generators of a gauge transformation can be 
constructed for the Hamiltonian using the so-called algebra of the Dirac constraints (l62ll 
(unnumbered equation preceding Eq. (29) of |lO[|)^^ 



r 1 



(127) 



These lead to the generator (see Eq. (29) of [lo| ) 

G = - j {[e^((H± + + {NTlie'') . + (nAT^) .) + e^H 



(128) 



+ 



m + ATf.H, + (AT^H,) . + N^.Ii) + ^H 



(129) 



There are some differences between these equations and the corresponding ones of [10| as 
we use Dirac's notation, e^**, which in ADM is called the three-dimensional metric. Actually, 
this renaming is sloppy and confusing because (see 01221] ) '^g^^ = ^ gkm _ _ s^fcm _|_ 

^ which, when solved for ^g'^"^^ is equivalent to Dirac's e^"^ f|T5|) . We keep e^'* to avoid 



From this point we are considering the ADM formulation and do not use the superscript ADM. 
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the temptation to raise spatial indices with this tensor. We can do this for the spatial metric 
(as they are inverses) but we cannot do this for derivatives, as was done in [l^ in the second 
term of fll29p where N'^Ylj is correct only if we consider A^'-' as short for N^ie^K This is 
because, according to the standard rules of raising indices in GR, 

= g^t^d, = g^'d, + = e^'d, + ^^d, + g^'^do, (130) 

and 



= e^'di (131) 

only if g^'^ = which is Dirac's simplifying assumption (1321) . If we use (11311) instead of (I130p 
we would obtain a different result. 

The generator (I129p allows one to find the transformations of the ADM fields (A^, A^* 
and gkm) and then by using the definition of these variables (11191 1120p , we can formally find 
the transformations of (7^1,. We do this in natural order - we first find the transformations 
of the ADM variables using the generator in terms of the ADM variables and then revert 
to the metric tensor. Transformations of the ADM fields are calculated using 6 {field) = 
{field,G}. 

Starting with the simplest variable, A^, we obtain 

SadmN = {N, G} = E^jN^ -ej,- s'N^i (132) 

— 1 /2 

which using A^ = {—g^^) gives 

SadmN =^-{-g''y'^'6ADMg'' 

and so we find 



SADMg'' = 2 i-g'Y'^' 6N = 2 {-g'^'^' 
This differs from the diffeomorphism transformation (11090 : 

X „00 _ r,fO,0 fO^OO ffc^OO 

0(diff)9 - 24 - 4 5',o - ? 9,k ■ 



[-9 '19 j - £,0 - ^ ( (-^ j 



no\-i/2 



133) 
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In ADM variables we cannot restore diffeomorphism invariance; the most that can be 
done is to present (11331) in a form similar to a diffeomorphism: 



,00 



and using = g^^d^j, = g^^do + g^^dk this becomes 



e- 



-L /' ^00\+l/2 



00\+l/2 



00 

y,o 



Ok 



g 



00 

g,k ■ 



The combinations in square brackets "correspond" to the diffeomorphism parameters 



(134) 



k _ , g^' ^± ( ^00\+l/2 ^^k _ ^ J 

N 



or 



(135) 



(136) 



lOO 



Equations fll34til37|) are equivalent to the result of 19|] (where the methods of 12j were 



n 



(137) 



used) and also to what is found in |2lL |2J]. The relations (11341 1136^ can be found in the 
Appendix of Castellani's article [l^, but (11351 11371) were not given there explicitly. This 
field-dependent redefinition of gauge parameters provides a "correspondence" 13], but not 
an equivalence with the diffeomorphism, that follows directly from consideration of the Dirac 
Hamiltonian. 

For the next variable, A^'^, we obtain 



SadmN'' = {A^^ G} = - U^Nje^'' - ef^Ne'"' + e^N^^ - e'^^N^ + e% 
and using A^'^ = — ^ir we find 



(138) 
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Ok 

gOO 



After expressing and A^'^ in terms of the metric 



^ADhig 



+g 



00 



(_^oo)-i/2j ^jk _ ^± ^_gOoyy\k, ^ i_g_ 



'J' ( ) + 



we again have an invariance that is not a diffeomorphism (llU9p 



0(diff)g =K +4 -K gfl ~K g,r. 

If we perform the field-dependent change of parameters fll36[ 11371) we again can present 
^ADhig^^ in the form of a diffeomorphism transformation. 

The transformation of the space-space components gkm was not considered in |10|] be- 
cause, according to the author, it is well-known that (11291) generates a diffeomorphism 
transformation of gkm- Let us check this statement; 



^gkm = {gkm, G} = [£^T-(-± + e'Ti-i] 

which, keeping only the fl'^'-dependent part of secondary constraints, gives 



(139) 



- {-n {g^v9^, - \g.,gv)j 2ALn^^ - 2 {e^g.,) ^ All + i^9p^,, - g,,.) Kl- 

(140) 

We must express n^"^ in terms of gtj and its derivatives; using Eq. (7-3. 9b) of [61] 



which is (taking into account that the "three-dimensional quantity" g^^ in ADM is Dirac's 

gip) 



rsabj^O 
ah' 



141) 
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This expression is equivalent to Dirac's expression for ( ITSll . Substituting (11411) into 
f fm and using ([HD, ([I9]) and ([20]), we obtain 



or in the exphcit form, using ( |T8l) . for the ADM formulation 



pq 
km'' 



^]m9ik — £]k9im " gkm,i (142) 



and for Dirac's variables 



^ADMSkra — {^9^^) ^ [d^^ {9a0,b + fi'60,a " Qabfl) + 9^^ i9ak,b + 9bk,a — 9ab,k 



~ ^]m9ik — ^]k9im — £^9km,i- (143) 

This is again different from the transformation of the spatial components of the metric under 
a diffeomorphism (11 101) 



S{diff)9km — -fi'fcm,0^° - 9k0^^m " 9m0^^k ~ 9km,iC " 9kiC,m ~ 9miC,k- (144) 

Again, only after the field-dependent redefinition of parameters (11361 1137p we can obtain 
a "correspondence" between SADM9km and 6(^diff)9km- (Note that both parameters have to 
be redefined despite the apparent equivalence of the last three terms in bot h eq uations (1143P 



and (I144p .) So, as in the case of Dirac's Hamiltonian, both methods 
same result (121) for the ADM Hamiltonian. 



Id, 



12 



produce the 



We would like to note that even a spatial diffeomorphism does not follow directly from 



the ADM formulation despite what is often stated in the literature (e.g. 29|]). If we treat 
the lapse and shift functions as 'multipliers'^^ and consider the secondary constraints as 



The authors of in the "Historical remark" on p. 486 stated that "The great payoff of this work [ADM] 
was recognition of the lapse and shift functions of equation (21.40) [the same as in (21.42) or our (|122p ] 
as Lagrange multipliers, the coefficients of which gave directly and simply Dirac's constraints." As we 
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being primary (the contradictions that result from such manipulations have already been 
discussed), we do not correct this problem because a spatial diffeomorphism does not follow. 
In this case, according to Castellani's procedure, the generator is simply 

G = e^H^ + e'Hi (145) 

which is equivalent to (11391) . Even in such a "formulation" the gauge parameters have to be 
redefined and, in addition, the transformations of lapse and shift functions equal zero (see, 
e.g. dm])). 

There is only one way to "derive" spatial diffeomorphism invariance and it explains the 
origin of the term "diffeomorphism constraint" . It behooves us to warn the reader that such 
a "derivation" has nothing to do with any procedure. If, in addition to eliminating primary 
constraints and promoting secondary to being primary (which leads to the generator (11451) ). 
we also consider only the second term of this generator, then 



^Qkm = {gkm,s'Hi} (146) 

will give the spatial diffeomorphism (see the second line of fll43p ). The only possible expla- 
nation of why such manipulations were accepted is that it seems to follow the "guidance" 
which comes from linearized gravity. From the derivation of the gauge transformations for 
linearized gravity in 371], it is clear that the only part of the generator proportional to x^"' 
contributes to the transformation of the space-space components of the metric tensor, but 
this is not the case for full GR. 

Ironically, the result (11461) . which is just the consequence of a series of manipulations that 
contradict any consistent procedure, is often presented as being the "problem" of the Hamil- 
tonian formulation of GR: "Hamiltonian" and "diffeomorphism" constraints are treated in 



a different manner 



29 



84| . Furthermore, the conclusion is drawn, based on (11461) . that 



"the diffeomorphism constraint can be shown to be associated with the invariance of general 
relativity under spatial diffeomorphism" 29|] (see also 7l|). Finally, because (I146P leads to 



have shown, "the great payoff' of this recognition is that the ADM formulation lost the connection with 
GR and the gauge transformations derived from it are different from diffeomorphism. We would also like 
to mention that in the same "remark" the authors wrote: "Dirac paid no particular attention to any 
variational principle". The interested reader is encouraged to look at Dirac's papers, especially at Q|, to 
recognize that this is not correct. 
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a spatial diffeomorphism, this result is interpreted as "disappearance of DiffA^" and con- 
sidered as "the problem that has worried many people working in geometrodynamics for so 
long" re- 
statements similar to the few of the forgoing quotations above can be found in many 
articles. They are based on questionable manipulations; but, at the same time, they clearly 
demonstrate that some authors correctly consider that the restoration of diffeomorphism 
invariance for all components of the metric tensor to be expected (however, "this expectation 
has never been fully realized..." [3])^^ and its absence is taken to be a deficiency or a 
contradiction arising in the Hamiltonian formulation. To have the possibility of restoring a 
transformation (whatever it may be) of all variables one must work in the full phase space 
and use all of the first-class constraints. We disregard "approaches" leading to 01451) and 
(11461) and return to the analysis of the total Hamiltonians (11151) and (I116p . 

We now ask why, in the ADM case, we must redefine the gauge parameters (I134til37]l to 
have a "correspondence" with diffeomorphism, whereas from the Dirac Hamiltonian (and 
also from the formulation without Dirac's modifications [s^) diffeomorphism arises directly. 
Dirac's Hamiltonian was obtained from the Lagrangian of GR, after some integrations that 
do not affect the equations of motion; and the ADM Hamiltonian apparently follows from 
the same Lagrangian. The two Hamiltonian formulations of the same Lagrangian ought to 
be equivalent and should give the same gauge transformation, which is not the case for the 
Dirac and ADM approaches. 

It is well-known that different sets of phase space variables can be used to describe 
iamiltonian systems. In the ordinary Classical Mechanics of non-singular systems (e.g., 
ItsI). for a given Hamiltonian H {qi,Pi) and the Hamilton equations, 

dH dH 
qi = {qi, H} = — , = {p„ H} = - — , (147) 
opi oqi 

we can pass to another set of phase space variables {Qi, Pi): 



qi = qi {Qk, Pk) , Pi=Pi {Qk, Pk) , K {Qk, Pk) = H {qi {Qk, Pk) , Pi {Qk, Pk)) , (148) 



such that 



This expectation was fuUy reahzcd in 30] and in the previous Section. 
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OK OK 
Q, = {Q„ K} = — P, = {P„ K} = -^^. (149) 

If this system of Hamilton equations can be solved, then one can return to the original 
variables by using the inverse of fll48p 

Qi = Qi {qk, Pk) , Pi = Pi {qk, Pk) ■ (150) 

In Classical Mechanics, even for non-singular systems, it is well-known, that [tJ] "...at 
first sight we might think that arbitrary point transformations of the phase space are now 
at our disposal. This would mean that the 2n coordinates q^ and Pi can be transformed into 
some new Qi and Pi by any functional relations we please. This, however, is not the case." 
For the non-singular systems, the necessary and sufficient condition that the transformation 
( 1148P to the new set of variables {Qi, Pi) is a. canonical transformation (and keeps the two 
formulation equivalent), is 

{Qi, Qkjg^p = {Qi (p, q) , Qk (p, q)}p^g = o, 

{P., P,}q^p = {P. (p, q) , Pfc (p, g)}^^, = 0, (151) 



{Qi, Pk}Q p = {Qi (p, q) , Pk {p, q)} 



PA 



Jik- 



The change of phase space variables (canonical transformations) for unconstrained Hamil- 
tonians is an old and well established topic that can be found in many textbooks on Classical 
Mechanics. For constrained Hamiltonians the situation is different and even the number of 



papers (e.g. see [76|, l77|]) that discuss the general questions or particular examples (which 
are mainly final dimensional and artificial) of such changes is minuscule compared to the 
number of articles in which such changes are used without any analysis of their consequences. 
Such changes are especially common in Hamiltonian formulations of OR (in both Einstein 



6| and Einstein-Cartan |78l. l79l. l80f| forms). 

The Dirac's generalization of the Hamiltonian formulation to constrained systems leads 
to the system of equations which is similar to (11471) : 

qi = {qi, Ht} , Pi = {pi, Ht] , 
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where Ht also includes the primary constraints. We restrict our discussion to gauge invariant 
systems with only first-class constraints. For such systems, the conditions that should be 
made on possible changes of variables in phase space seems to be more restrictive. The reason 
for this is that gauge invariance depends on all of the first-class constraints and their PB 



algebra 



10 



11 



12l |. In addition, it is related to the total Hamiltonian of a particular model. 



in contrast to unconstrained systems where the conditions (11511) are, in fact, independent of 
the Hamiltonian. For gauge invariant systems the change of variables fllSOp must preserve 
gauge invariance, i.e. gauge invariance derived in the {Qi, Pi) variables, after using the 
inverse transformations, must produce the same result as in the original {qi,Pi) variables. 
The set of phase space transformations fll48p that preserves this property is the equivalent 
set or, as in non-singular case, we can call such changes canonical transformations. 

In linearized versions of the Dirac formulation (see the previous Sections) and that of 41 
which were considered in [37[, both formulations, despite having different constraints and 
Hamiltonians, lead to the same gauge invariance. The relation between the two formulations 
was discussed and it was shown that they are related by canonical transformations (11511) . 
which is exactly the condition known to be needed for non-singular Hamiltonians. Despite 
there being far more complicated expressions for constraints and transformations, a similar 
relation (11511) exists between the non- linearized formulation of Dirac and SOj; and these, as 
we have shown, also have the same gauge invariance. These examples demonstrate that the 
ordinary condition for the transformation to be canonical, which is known for unconstrained 
Hamiltonians is also correct for the general (constrained) case; i.e. it is a necessary condition 
as before. Yet, it is not sufficient and we shall demonstrate this fact by way of example. 

Surprisingly, we have been unable to find, in the literature, either a discussion of the 
relations between the ADM and Dirac phase space variables or the transformations between 
them; which is strange as many authors presume their equivalence by calling this formulation 
"Dirac- ADM". In 81| the authors called the variables A^, A^* and gkm "an equivalent set" 
which "is analytically convenient and geometrically more significant". In footnote 5 of 81 1 
(appeared in 1959) it is stated that "The properties of these variables are discussed in detail 
in a forthcoming paper by C. Misner" that we have not been able to find. The convenience 
and significance are not our present concern, but the question of equivalence of the ADM 
and Dirac sets of variables is important. We are interested in the equivalence of gauge 
transformations in the two approaches, i.e. we must work in the full phase space. The 
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complete relation between these two sets of variables is not known, but there is one PB, 
(x) , n'^' (a;') }, that can be easily checked. The space-space components of the metric 
tensor and corresponding momenta are the same in both formulations, i.e. H*^' = II'^' (p'^') = 
p^K The n^' of ADM given by (11411) is equivalent to p^'' of Dirac given by flTH]) . It is sufficient 
to check using the Dirac variables g^^, p'^^ the PB (x) , II'^' (a^')} which, if ADM variables 
are canonical, must give zero. Using the corresponding fundamental PBs (11 24111 26|) we 
obtain: 



} 9,P OQkl ^ 

(152) 

Once again, we have the result that depends on Dirac's simplifying assumption ( l32i) : if 
we impose g^^ = 0, then the PB of (11521) gives zero. In general this PB is not zero and the 
transformation from {g^uiV^") (^^n), (A^^IIj), and {^gkmi^^^^ is not canonical. One PB 
is enough to show this, irrespective of the results we might obtain for the PBs among the 
other phase space variables. 

We take a note that equation (I152p gives zero in the static coordinate system, but in this 
case the corresponding components, gok, and their conjugate momenta have to be dropped 
out of the formalism from the beginning, as in the case of the Hamiltonian formulation in 



the Schwarzschild metric 



391] where only four components of the metric tensor are left in 



the Lagrangian before passing to the Hamiltonian. 

This simple calculation, ( ]152p . allows us to conclude that the two Hamiltonians of Dirac 
and ADM are not related by a canonical transformation and the respective failure of the 
ADM variables and the Hamiltonian to produce a diffeomorphism transformation is a mani- 
festation of this non-equivalence. Moreover, (11521) shows that the ADM variables are not the 
canonical variables of GR. The converse statement is also true and the metric tensor is not 
a canonical variable of the ADM formulation. The ADM formulation might be considered 
as a model (geometrodynamics or ADM gravity) without any reference to Einstein GR, 
but in this case a "correspondence" between the two transformations (1134111371) is, in fact, 
meaningless. The transformations that follow from ADM are given by (I132p . (I138p . and 
(11421) and in the absence of canonicity we cannot return to the transformations of the metric 
tensor. There is another characteristic that supports the loss of connection with the original 
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variables: it is impossible to find any redefinition of 11 and 11^ in terms of Dirac's phase 
space variables (whether they satisfy fllSip or not) that can transform his total Hamiltonian 
(dUD into the ADM total Hamiltonian (ITT^ . 

In general, no algorithm exists for finding a canonical transformation but the canonic- 
ity of a given transformation can be checked. There are no canonical transformations for 
Hamiltonians that involve only a change of generalized coordinates, as this change must be 
accompanied by transformations of the momenta that can be found by using the following 
procedure [t^. If the transformations of the generalized coordinates (fields) are given, one 
can find the corresponding transformations of the momenta that will guarantee that the new 
coordinates and momenta are canonical and satisfy (11511) using the relation 

PiSqi = PiSQi. 

Let us see what we can obtain from this relation for the ADM change of variables 
By performing the variations 6Q = g^R using (I119til20p . we find that 



which gives 



(^9af3 dg^p dgafs 



p'' = -T^l{-9n'^\ (153) 



p""- = ni i-g'')-'^' g'- + U,^e'- (154) 



and 



(155) 



Now solving for the II's: 

n = -2(V°)"'^'/°, (156) 
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n„ = 2g,„„p°'» + 2g(k,p»°, (157) 



n"=P" + ?5T^P"-!BP°'-!5TP°'' (158) 



Oq Op „0p „0q 

I y 00 _ y Oq _ y_ 

gOO gOoP gOoP gC 

Note that, to have 11^^ = j/"', as in the case of the Dirac and ADM Hamiltonians, we again 
must impose the condition g^'' = 0. 

Only if equations fll56til58l) are taken together with the relations for the generalized 
coordinates (111911120]) . are the transformations canonical. It is not difficult to check that the 
canonical properties of PBs are preserved, and as an example, the PB that we considered in 
ffT^ gives 

{N (.) ,n« (.')} = { i-rr" y + p^v"" - ^p°' - p°']^= °. 

as it should for variables that are connected by canonical transformations. For non-singular 
Lagrangians and their corresponding Hamiltonians, the ADM change of variables accom- 
panied by the change of momenta (1153111551) would be sufficient to obtain the new set of 
canonical variables. If GR were a non-singular theory, these transformations would guaran- 
tee equivalence between the two formulations; but for the constrained Hamiltonian this is 
not the case. If the canonical transformation of (11221) . (I153H155|) are performed in the Dirac 
Hamiltonian, we will not obtain the ADM Hamiltonian, and we will not obtain a consistent 
result. In particular, for a description of a constrained system, the total Hamiltonian, Ht, 
is important as it includes all the primary constraints. In the Dirac formulation, these are 
(I22D: 

^7oo,o/° + 2(7ofc,o/'. (159) 
Substitution of (11221) . (1153111551) into this equation gives 

iV,oH + N'^U^ + g.j^o ( 2^^^ + Nm^e^^ 

which is nonsensical. The first two terms are equivalent to (11591) but an extra term appears, 
which is zero only if A^'^ = 0. The space-space velocities have already been eliminated in 
favour of their corresponding momenta but now they reappear and it is not clear what to 
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do with them at this stage. If we treat them on the same footing as the rest of the terms 
with time derivatives, we must specify their coefficients as primary constraints, which would 
then give a total of ten primary constraints. The same change of variables in the canonical 
part of the Hamiltonian will give contributions that are quadratic in all momenta. Without 
further analysis we see that this will probably lead to some contradictions, as it is clear that 
the constraint structure of the Hamiltonian is changed and one would expect second-class 
constraints, etc. 

The foregoing example, flll9H120p and (1156111580 . clearly demonstrates that the condition 
(11511) . which is necessary and sufficient for the transformation to be canonical in the case of 
non-singular Lagrangians, is not a sufficient for singular Lagrangians. 

Actually, for the ADM change of variables, the non-canonical nature of the transfor- 
mations is immediately clear and this is not even related to a singular structure of the 
Lagrangian of GR. The possible existence of additional restrictions beyond (11511) is under- 
standable, and for the Hamiltonian with first-class constraints it can be expected. Gauge 
transformations are derived from the first-class constraints and the whole PB algebra of 
constraints plays a key role in this derivation. It is not enough to have the same number of 
constraints in the two formulations. The ADM change of variables keeps the same number 
of constraints as the Dirac formulation; but the PB algebra of constraints is affected. The 
simplest example is a PB among primary and secondary constraints, which are zero in the 
ADM case and proportional to true constraints in the case of the Dirac Hamiltonian (|63ll . 
From a mathematical point of view (we have already spent enough time on the interpreta- 
tional aspects) the ADM formulation is just the result of a non-canonical change of variables 
in Dirac's Hamiltonian of GR and if Dirac's formulation allows one to derive the diffeomor- 
phism transformations, then the ADM formulation, because of this non-canonical change of 
variables, does not allow one to restore either the full diffeomorphism invariance for all com- 
ponents of Qfj^u or even for its spatial part (as it is usually claimed) without a non-covariant 
and field-dependent redefinition of gauge parameters. In addition, we observed that some 
equations in the ADM formulation are true only if gok = 0. So the ADM change of variables 
is somehow also related to a static coordinate system, but in a strange and obscure way: 
gok is not zero at the outset, but later should be set equal zero so that some subsequent 
relations are made valid. 

We conclude that the Dirac Hamiltonian of GR was obtained by following the "rule 
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of procedure" and, because of this, it is canonical already and the use of the adjective 
"canonical" is a tautology. It is not equivalent to the ADM Hamiltonian, which, as we have 
demonstrated, is the result of a non-canonical change of variables. The ADM formulation is 
obtained by abandoning the "rule of procedure" and, consequently, only by a canonization 
can it be ironically called the "canonical formulation of GR" . 

One can argue that the ADM Hamiltonian is not obtained from the Dirac Hamiltonian. 
There are no references in jo] to Dirac's article [sl, the only reference where a derivation of 
Dirac's Hamiltonian has been considered. In their culminating paper jo] and in 13 preceding 
articles it is mentioned only once in 72| and in a different context. However, according to 
both Dirac and ADM, their respective Hamiltonians are derived from the same theory - GR. 
In both cases some modifications of the EH action were performed. In Dirac's case these 
modifications are explicitly stated, in the ADM case it is more difficult to trace what has been 
done. Let us start from the Dirac Lagrangian. We have discussed above how, in the course 
of a Hamiltonian analysis that follows Dirac, the possibility of additional integrations by 
parts appear when the Hamiltonian is to be expressed as a linear combination of secondary 
constraints (see our discussion after ( l55i) ). The modified Lagrangian can be written in the 
following form (up to total derivatives as in fl5^ ) 



LDirac = ^ a/ det gkm ("5'°°) E'' 



-,Ofc 



~ 9ab,0 H i9ak,b + 9bk,a — 9ab,k 



X 



9l\ _ ^9^( 

9sq J 9rs,0 1 \9rm,s + 9sm,r 9rs 



+ Vdet gkm (-^°°) 



00\-V2 



^7m„,Hi5""'* + \9mn,k9p,,t [E^'^^'^ 6^' - 2E'=*P"e"^« - AE^'^'^' 6^^) 



(160) 



By performing the change of variables of fll22p in (I160p . we obtain the ADM Lagrangian 



^ADM 



7V^ 



9kr, 



N 



:E 



rsab 



{9apNP)^b + i9bpNP)^a ~ 9ab,0 - N'' {gak,b + 9bk,a " gab,k) 



{grqN'^) s + {9sqN'^) r " 9rs,0 " N"" {grm,s + gsrn,r " 5'rs,m) 
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+ a/ det gkmN 



(161) 



Or, by using the intrinsic and extrinsic curvatures, '^Rrs and Krs-, respectively, it can be 
written in more familiar form HQ 



Ladm = V^^Q^N (^"'^^iT^.i^,, +^ R) (162) 

with 

Krs = {Nr\s + Ns\r - 5'rs,o) 

where, as before, "|" means covariant derivative with respect to three dimensional metric 

{Nrls = Nr,s-T';,Nk). 

Now using fll6ip or (11621) we can easily obtain the ADM Hamiltonian, which is the same 
as (11151) . However, following such a detour we cannot avoid the question of whether the 
ADM variables are canonical. The transformations of the metric tensor derived from the 
Hamiltonian formulation of the Dirac Lagrangian (11601) and from the Hamiltonian formula- 
tion of the ADM Lagrangian (11611) are different, the destination is changed. This detour is 
a wrong turn or perhaps a dead end road. 

We combine the results of these two formulations into a compact visual form^^: 

-'-'Dirac 
H Dirac 

{p,q) Hadm{P,Q) 
^ ^ (163) 

SadmQ 

SdiffQ ^ SadmQ 



It is reasonable to expect that if by a change of variables =^ we can obtain a new equiv- 
alent Lagrangian and find the corresponding Hamiltonian, then the two Hamiltonians (in 

This 'pictorial visualization' is based on results of calculations, not the other way around. 
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both the new and the original variables) should also be equivalent, be related to each other 
by a canonical transformation, and necessarily lead to the same gauge invariance. If the 
Hamiltonians are not related by a canonical transformation, then the two Lagrangians are 
not equivalent. This is the natural conclusion that one can make. In particular, application 
of the Lagrangian methods used by Samanta [26:] to derive the diffeomorphism invariance 
of GR, when applied to the ADM Lagrangian, should not give the diffeomorphism transfor- 
mations (82| but rather the same transformations as obtained in its Hamiltonian treatment 
in 191]. This is a consequence of the well-known equivalence of the Lagr ang ian and Hamil- 



17 



: the "vertical" 



tonian formulations for any system, either non-singular or singular 
equivalence of (11631) . If a "horizontal" equivalence is broken, either for the two Hamiltonians 
or for the gauge transformations, then it is broken everywhere, including at the Lagrangian 
level. 

The Hamiltonian formulations of the linearized versions of the Dirac and gamma-gamma 
Lagrangians lead to the same algebra of constraints and gauge transformations (although the 
Hamiltonians themselves and the constraints are different). But the two Hamiltonians are 
related by a canonical transformation 37|. Such transformations also exist in the case of the 

the PB algebra of constraints can be easily seen by comparing [30| to the results of Section 
II. The explicit canonical transformation and its effect on constraints, their algebra and 
structure functions will be given in 42]. 

The main subject of this article is the Hamiltonian formulation of GR in second-order 
form and, in particular, a comparison of the formulations related by a change of variables. 
However, we think that some comments on similar changes made to its Lagrangian should 
be given. 

It is often stated that in a Lagrangian formulation of a model any field redefinition is le- 
gitimate provided it is invertible (i.e. it has a non-zero Jacobian). For singular Lagrangians, 
and especially gauge invariant ones, this is obviously not a sufficient condition. One addi- 
tional restriction in gauge invariant cases is the preservation of the rank of Hessian as this 
gives us the number of gauge parameters (if all constraints are first-class). We cannot have 
two equivalent formulations if they have a different number of gauge parameters, and we 
cannot, by a change of variables, eliminate some gauge invariance or create a new gauge 
invariance. If we were to make such a change, we can of course, treat the new Lagrangian as 
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some different model, but we cannot relate it to the original one as any connection with the 
original theory is lost. Obviously these two conditions are necessary, but not sufficient. The 
ADM charige of variables satisfies them both, but leads to different gauge transformations 
(compare SOj and Section III versus [19] and (1134111371) ). A change of variables in singu- 
lar (in particular, gauge invariant) Lagrangians is a much more restrictive procedure if one 
intends to preserve its equivalence with the initial formulation. One way is to rely on the 
Hamiltonian method and check if the new variables are canonical and the two total Hamil- 
tonians are equivalent, including the primary constraints. We must also check whether the 
entire algebra of constraints is equivalent, as this algebra is responsible for the gauge trans- 
formations. This can be considered as a confirmation of Dirac's statement 7|] "I feel that 
there will always be something missing from them [non-Hamiltonian methods] which we can 
only get by working from a Hamiltonian, or maybe from some generalization of the concept 
of a Hamiltonian". However, we think that some criteria for the equivalence between two 
sets of variables for singular Lagrangians can be formulated at the pure Lagrangian level. 
At the Lagrangian level, a gauge invariance is related to the existence of gauge identities 
and an inappropriate change of fields can modify or even destroy them. This echoes the 



conclusion of Isham and Kuchar 73| that "... space-time diffeomorphism has somehow got 
lost in making the transition from the Hilbert action to the Dirac-ADM action" 

There is another indication of the incorrectness of the ADM change of variables at the 
Lagrangian level that comes from Numerical Relativity. In almost all methods of numerical 
integration of the Einstein equations the starting point is the ADM 3+1 decomposition. So 
at the outset the Einstein equations are replaced by the ADM equations [6]. It was shown 



that, in contrast to the Einstein equations which are strongly (strictly) hyperbolic (SH) |85| . 



86 
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This change in the 



the ADM equations are weakly hyperbolic (WH) (e.g., see 
type of equations is related to the different constraint structure and different transformations 
derived from a non-equivalent Hamiltonian. From a computational point of view there is a 
fundamental difference between SH and WH systems of PDEs: the former are well-posed 
and convergent, whereas the latter are not well-posed and divergent 89|. There is numerical 



evidence that ADM-based algorithms are unstable. As is indicated in 



90|: "The common 



We have demonstrated the inequivalence of the Dirac and ADM formulations. In 73| the authors discussed 
geometrodynamics which uses ADM variables and their statements should be applied to the ADM action 
only. 
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lore these days is, however, that the standard Arnowitt-Deser-Misner (ADM) formulation is 



the one which most easily suffers instabilities". Or in [89|]: "It took several years to realize 
that such instabilities were not associated with the numerical algorithms but rather with 
the mathematical structure of the ADMY [ADM [6] and York 9l|] system itself. This is 
not just an additional indication of the incorrectness of the ADM change of variables, but 
also a demonstration that even as a model, geometrodynamics is an ill defined theory. 

In general, the change in the type of equations, from SH to WH, or a change of "level 
of hyperbolicity" 9^, is an indication that in the process of transforming from the Einstein 
to the ADM equations, some "damage" was done (as there is no longer a complete set 
of eigenvectors associated with the characteristic matrix 93]). The proposed "cure" 9^ 
of such "damage" in most approaches lies in a modification (or "adjusting") of the ADM 
equations by adding terms involving constraints [9J] (trying to restore what has been lost) 
or by using different choices of the lapse and shift functions j90| to give the ADM system 
of equations well-posedness (or "quasi well-posedness" (90|]), rather than returning to the 
original Einstein equations and constraints that preserve diffeomorphism invariance. 

The main point of this relatively long Section is not to prove that the ADM variables are 
not canonical variables for GR (as is shown by the one simple PB fll52p ). but to demon- 
strate and discuss the restrictive conditions that must be made on change of variables in 
any Hamiltonian formulation of a singular Lagrangian, using GR as an example. A blind 
change of variables in singular systems without performing a thorough analysis and without 
developing mathematical criteria for such changes can lead to a wrong result. All new vari- 
ables that are introduced in such cases, regardless of their physical or geometrical meaning, 
regardless of what new names were given to reflect their interpretation, or after whom new 
variables were named, must be carefully analyzed if one wants to keep all the properties of 
the original theory intact or, in other words, if one intends to study the original theory and 
not a substitute, which can be ill defined, even as unrelated to the original theory model. 
One should under no circumstances attribute to the original theory the contradictions or 
problems that arise after such inappropriate changes are made; and never project any novel 
result or discovery obtained by abandoning a "regular and uniform rule of procedure" into 
the original theory or to Nature Herself. 
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